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Abstract 

The Skyrme model is reconsidered from an effective theory point of view. From the most general 
chiral Lagrangian up to including terms of order p , N c and 5m 2 (6m = m s — m), new interactions, 
which have never been considered before, appear upon collective coordinate quantization. We 
obtain the parameter set best fitted to the observed low-lying baryon masses, by performing the 
second order perturbative calculations with respect to 8m. We calculate the masses and the decay 
widths of the other members of (mainly) anti-decuplet pentaquark states. The formula for the 
decay widths is reconsidered and its baryon mass dependence is clarified. 
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I. INTRODUCTION 



Evidence of a new baryonic resonance state, called + (154O), has been claimed recently 
by Nakano et. ai.Q, with S — +1 and a very narrow width T < 15 MeV. Several other 
experimental groups have confirmed the existence P, 0,0]. It appears in the recent version of 
Reviews of Particle Physics^ with the *** rating, though itsparity has not been established. 
Evidences of less certain exotic pentaquark states, 0(1860) p and O^" [3, Q] have also been 
claimed. The discovery of the pentaquarks is expected to lead us to a deeper understanding 
of strong interactions at low energies. In reality, it stimulates new ideas and reconsideration 
of old theories and experimental data. 

The discovery was motivated by a paper by Diakonov, Petrov and Polyakovj^J. They 
predicted the masses and the widths of the anti-decuplet, of which + is presumed to be a 
member, within the framework of the "chiral quark-soliton model" (xQSM) 1 . See Ref. j^| for 



a review of the xQSM. The chiral quark-soliton model prediction is reexamined in Refs. 



See also Refs. 



Jaffe and Wilczek 



16 



113, 



~y\ for earlier papers on the anti-decuplet in the Skyrme model. 
18| | , on the other hand, proposed a quark model picture of pentaquarks 
based on diquark correlation. See Ref. 3] for a review of this approach and other interesting 
aspects. 



Since Witten pointed out that baryons may be considered as solitons|25j] in the large- N ( 



limit 



n 



26], and showed that the soliton ("Skyrmion") has the right spin-statistics |27[ thanks 

to the Wess-Zumino- Witten (WZW) termHQ, 

much effort has been done to explore the 

consequences. But the results do not look very successful. In order to fit the results to 
the observed values of masses, coupling constants, such as the pion decay constant, which 
appear in the chiral Lagrangian must be very different from the experimental values. For 
example, in the SU(3) Skyrme model, the pion decay constant becomes typically one third 
of the experimental value to reproduce the correct mass splitting bOjl. It also predicted an 

nnn 

anti-decuplet |15L I16L I17|. which many did not believe to exist at that time. 

But + has been discovered! It is time to take a serious look at the Skyrme model again. 
One of the most important aspects of G + is its narrowness. Several analyses of older data 



1 The xQSM has its own scenario based on instantons. For our purpose, however, it is useful to regard it 
as a version of the Skyrme model jicj with specific symmetry breaking interactions. See for a recent 
review on the background of the xQSM. 
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Interestingly, the Skyrme 



indicate that the width may be less than 1 MeV[20j, |21 
model is believed to be capable to explain the narrowness. It is claimed that the width 
becomes very narrow due to a strong cancellation [2^. 

A natural question is: Is this a general result of the Skyrme model, or a "model- 
dependent" one? What is the most general Skyrme model? The Skyrme model is noth- 
ing but a model. But, as Witten emphasized, the soliton picture of baryons is a general 
consequence of the large- N c limit of QCD. Assuming that the large- N c QCD bears a close 
resemblance to the real QCD, we may consider an effective theory (not just a model) of 
baryons based on the soliton picture, which may be called as the "Skyrme- Witten large- iV c 
effective theory." 

Our key observation is that the Skyrme model conventionally starts with the particular 
chiral Lagrangian, which consists of the kinetic term, the Skyrme term (which stabilizes the 
soliton), the WZW term and the leading SU(3) breaking term, 

S Sky r m e = d A xTr (d^Ud^) + / ^Tr ( [^U, tfd v U] 2 ) + N C T[U] 



FIB 



dVTr {M ] U + MU ] ) 



where Ai is the quark mass matrix 2 



M 



m 
m 
m, 



;i.2) 



and 



T[U] 



240vr 2 



dZ ijklm Tr 



Q 



{diU) f/ f {djU) (d k U) {diU) f/ f (d m U) 



1.3) 



is the WZW term. The chiral perturbation theory |3lLl32j (%PT) is however an effective field 
theory with infinitely many operators 3 . We should keep it in mind that there are (infinitely) 
many other terms and the expansion must be systematic. 

In this paper, we explore such an "effective theory" approach, i.e., an approach based 
on a systematic expansion of the operators based on the symmetry and power-counting. 



2 We do not consider the isospin br eaki ng in this paper. 

3 Long time ago. Kindo and Yukawa[33( considered the Skyrme model in the framework of the xPT context, 
but their work did not seem to attract much attention at that time. 
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The parameters are determined by fitting the results to the experimental values. Once the 
parameters are fixed, the rest are the predictions. As an application of our approach, we 
give the predictions to the masses and the decay widths of the other (mainly) anti-decuplet 
baryons. 



Generalizations o 
ing vector mesons \34 



:he Skyrme model have been considered by several authors, by includ- 



37 



3811 . Inclusion of heavier vector mesons 



35l . |36| or the radial modes 
may correspond to the inclusion of higher order terms in the usual pseudoscalar Lagrangian, 
while the inclusion of the radial modes is beyond the scope of the present paper. 

The paper is organized as follows: In Sec. |HJ we derive the collective coordinate quan- 
tized Hamiltonian based on the effective theory approach. The Hamiltonian contains several 
SU (3) flavor symmetry breaking interactions which have never been considered in the liter- 
ature. In Sec. IIII1 we consider the eigenstates of the Hamiltonian, with the eigenvalues being 
the baryon masses. If the symmetry breaking terms were absent, the eigenstates form the 
flavor SU(3) representations. The symmetry breaking terms mix the representations. We 
calculate the masses up to the second order in the perturbation theory with respect to the 
symmetry breaking parameter Sm = m s —m. The masses are represented as functions of the 
parameters of the theory. In Sec. IIVI we numerically fit the calculated masses to the exper- 
imental values to determine the parameters. After determining the parameters, we obtain 
the masses of mainly anti-decuplet baryons, N', an excited nucleon with I{J P ) = |(| + ) ; an d 
E' with The decay widths are calculated in Sec. El after reconsidering the derivation 

of the formula for the widths in the collective coordinate quantization. Unfortunately, our 
calculation of the widths suffers from large ambiguity. We summarize our results and dis- 
cuss some issues in Sec. IVII The notations, conventions, and the derivations of several useful 
mathematical formulae are delegated to Appendix El Various matrix elements used in the 
calculations are summarized in Appendix [Bl A lot of tables are given there. The reason 
why we present them is that most of them have never appeared in the literature and it 
requires much labor to calculate them. The "traditional" approach to the Skyrme model is 
reconsidered from the new perspective in Appendix O The results obtained in Appendix O 
may be viewed as an evidence that our basic strategy is right. Finally in Appendix El we 
perform a parallel analysis with the xQSM symmetry breaking terms, 

H? pp = aDjg ^=D^F\ (1.4) 
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without considering how these interactions are derived. The parameters a, (3, and 7 are 
determined in a similar way, and the decay widths are calculated too. 



II. COLLECTIVE HAMILTONIAN 

A. Chiral Lagrangian up to including 0(p A ) and 0(N C ) 

Effective field theories are not just models. They represent very general principles such 
as analyticity, unitarity, cluster decomposition of quantum field theory and the symmetries 



of the systems |39p. The chiral perturbation theory (xPT), for example, represents the low- 
energy behavior of QCD (at least) in the meson sector. 

Although baryons in the large- N c limit behave like solitons, it is not very clear in what 
theory they appear. A natural candidate is the x?T, because, as emphasized above, it is a 
very general framework in which the low-energy QCD is represented. It seems that if baryons 
may appear as solitons, they should appear in the xPT, with infinitely many operators. 

At low-energies, only a few operators are important in the x?T Lagrangian. We may 
systematically expand the results with respect to the typical energy /momentum scale, p. 
This is the usual power counting in the xPT, and we assume it is the case even in the 
soliton (i.e., baryon) sector. 

To summarize, a general Skyrme-Witten soliton theory may be a systematic expansion 
of the soliton sector of the x?T, with respect to N c and p. Therefore, our starting point is 
the 577(3) ;\PT action (without external gauge fields) up to including 0(p A ), 

S xPT = S + S 1 + O(p 6 ), (2.1) 

5 = S J d A xTr (d^Ud^) + J d 4 xTi (M j U + MU ] ) + N C T[U), (2.2) 

51 = L x J d 4 x [Tr (dpUdPrf)] 2 + L 2 J d 4 xTr (d^d u U) Tr (d»U ] d v U) 

+ L 3 J d A xTr (dJJ ] d»Ud v U ] d v U) + L A B J d 4 xTr (d^d^U) Tr (M*U + MU ] ) 

+ L 5 B J d 4 xTr [dprtdPU (M j U + U ] M)] + L 6 5 2 J d 4 x [Tr (M j U + MU ] )] 2 

+ L 7 Bl J d A x [Tr (M j U - MU ] )f + L 8 5 2 J d 4 xTi (M ] UM ] U + MWMW) , 

(2.3) 



where Li (i — 1, • • • , 8) are dimensionless constants. The definition of these parameters is 



the same as that of Ref. 



except for Fq. Our normalization of Fn is more popular in the 



Skyrme model literature. The N c dependence of them is known |32t l40j]. 

B , 2Li - L 2 , L 4 , L 6 , L 7 • • • 0(iV c o ), (2.4) 
F 2 ,L 2 ,L 3 ,L 5 ,L 8 ••• ^(iV, 1 ). (2.5) 

In the following, we keep only the operators whose coefficients are of order N c . Experimen- 
tally, these constants are not very accurately known. We further assume that the constants 
Li, L 2 and L 3 have the ratio, 

Li : L 2 : L 3 = 1 : 2 : -6, (2.6) 

which is consistent with the experimental values, L\ = 0.4 ± 0.3, 1L\ — L 2 = —0.6 ± 0.5, 
and L3 = —3.5 ± 1.1 (times 1O _3 )0]. Note that vector meson dominance also implies this 
ratio [42]. This assumption simplifies the analysis greatly, due to the identity, 

Tr (ABAB) = -2Tr (A 2 B 2 ) + ^Tr (A 2 ) Tr (B 2 ) + (Tr {AB)f , (2.7) 

which holds for any 3x3 traceless matrices A and B. By using it, the L 1; L 2 , and L 3 terms 
are made up to a single expression, 



-L- 2 Tr([U^,U,U^d u U} 2 ), (2.8) 

where we introduced L 2 = l/(16e 2 ). This term is nothing but the Skyrme term 4 . (If we 
would not assume these exact ratios among Li, L 2 and L 3 , we would have extra terms which 
lead to the terms quartic in time derivatives of the collective coordinates, and would make 
the quantization a bit harder. Because we consider the case in which the "rotation" is slow 
enough, such terms could be ignored.) 
We thus end up with the action, 

S[U] = S J d 4 xTi (dpUPrt) + ^ J d 4 xTi (AWd^tfdJjf) + N C T[U] 

J d 4 xTi {M ] U + MU ] ) + L 5 B J d A xTi (d^d^U (M j U + U ] M)) 
+ L 8 B 2 J d 4 xTr (M j UM j U + MU ] MU ] ) , (2.9) 

4 We do not know who first noticed this fact. Probably this is widely known. Wc learned it from Ref. [l^ . 



which is up to including 0(N C ) and 0(p 4 ) terms. 

It is important to note that, though (|2.9jl is our starting point, we always need to keep 
in mind that there are infinitely many higher order contributions. It is the symmetry and 
the power counting that actually matter. The discussion given in this section (and in 
Appendix [CJ) may be considered as a heuristic derivation, which is however very convenient 
because it explicitly shows the (leading) orders of N c , m and 8m for various parameters. 

In the usual %PT, one loop quantum effects of mesons can be incorporated consistently 
to this order and the physical parameters such as the pion decay constant F n are so defined 
as to include one-loop corrections. In our analysis, we ignore all the quantum effects of 
mesons, but there are still tree-level contributions to physical parameters from the higher 
order terms. For the decay constants, we have 

F w = F (l + (2m)K 6 ), (2.10) 



F v = F ( 1 + \{m + 2m s )K 6 ) , (2.12) 



where 



k = F (l + (m + m s )K 6 ), (2.11) 
2 
3' 

K, = ^L 5 . (2.13) 

Meson masses are obtained by looking at the quadratic terms in the Lagrangian when 
expanded around U = 1, 

Ml = B (2m) (1 + (2m)K 3 ) , (2.14) 

Ml = B {m + m s ) (1 + (m + m s )K 3 ) , (2.15) 

Ml = ^B (m + 2m s ) M + ~(m + 2m s )K 3 \ + K 5 , (2.16) 



where 



32.00 



K 3 = ^(2L 8 -L 5 ), (2.17) 

-Tn 



,512Bg 





K 5 = ( ms - m y -£L 8 . (2.18) 



B. Collective coordinate quantization 



In this subsection, we derive the Hamiltonian which describes the baryons by using the 
collective coordinate quantization. In this treatment, we consider the soliton as a "rigid rota- 
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tor" and do not consider the breathing degrees of freedom, though such "radial" excitations 
should be important if we consider other states, such as those with negative parity 5 . 

There are two important criticisms against the above mentioned treatment. The first is 
an old one that the flavor SU (3) symmetry breaking is so large that the perturbation theory 
does not work. The so-called "bound-state" approach has been advocated by Callan and 



Klebanov 



43jj and it was recently reconsidered after the discovery of the + resonance 



It is Yabu and Ando|45j who showed that the "exact" treatment of the symmetry breaking 
term gives good results even though the collective coordinate quantization is employed. 
Later, it was shown that the perturbation theory is capable to reproduce the qualitati vely 
equivalent results if one includes mixings with an enough number of representations jlTl. . 
The second is claimed by Cohen j^,^] that, in the large- N c limit, the "rotation" is not slow 
enough for the collective treatment to be justified. Diakonov and Petrovj^J emphasized that 
due to the WZW term, the "rotation" is slow enough even in the large- iV c limit. 

The collective coordinate quantization with the flavor SU(3) symmetry is different from 



22l 



the one with the isospin due to the existence of the WZW term. See Refs. 

0,0, HQ- 

The rotational collective coordinates are introduced as a time-dependent S , [/(3)-valued 
variable A(t) through 

U(t, x) = A{t)U c (x)A\t), (2.19) 
where U c (x) is the classical hedgehog soliton ansatz, 



U c (x) 



exp {it ■ xF{r)) 





(2.20) 



V 00 V 

with the baryon number (topological charge) B = 1. The profile function F(r) satisfies the 
boundary conditions F(0) = n and F(oo) = 0. By substituting Eq. (j2.19J) into Eq. ()2.9|) . 
one obtains the Lagrangian 

L = -M d + )<o a I a ^A)^ + ^=u 8 - V(A), (2.21) 



5 We assume that the pentaquark states have positive parity. Otherwise they would not be "rotational" 
modes of the soliton, and our analysis in this paper would not make sense at all. 



s 



where we have introduced the "angular velocity" uj a (t) (a = 1, • • • , 8) by 



A\t)A(t) = -J2^ a (t), (2-22) 

a=l 

with A a (a = 1, • • ■ ,8) being the usual Gell-Mann matrices. 

The first term represents the rest mass energy of the classical soliton. In the "traditional" 
approach, the ^-independent part is a functional of the profile function F(r). The classical 
mass M c i is given by minimizing (the minus of) it by varying F(r) subject to the boundary 
condition. One might think that the parameters in the xPT action (J2.9)) have already 
been given and M c i would be determined completely in terms of these parameters. It is 
however wrong because the action contains infinitely many terms and therefore there are 
infinitely many contributions from higher orders. We do not know all of those higher order 
couplings and thus practically we cannot calculate M c i at all. A more physical procedure is 
to fit it to the experimental value. This is our basic strategy in the effective theory approach. 
The operators are determined by the action, reflecting the fundamental principles and 
symmetries of QCD, while the coefficients are fitted to the experimental values. Because we 
do not know the higher order contributions, the number of parameters is different from that 
of ()2.9|) . For the comparison with the "traditional" approach, see Appendix which also 
serves as a "derivation" of the terms discussed below. 

The most important feature of the Lagrangian (j2.21jl is that the "inertia tensor" I a p(A) 
depends on A 6 . From the symmetry of the ansatz and the structure of the symmetry 
breaking, it has the following form, 

I aP {A) = I% + I' ap {A), (2.23) 

with 



1 af3 



Ix6 aj 3 a, (3 el 

hS a p a,(3ej (2.24) 
otherwise 



6 A mechanical analogy is a top with a thick axis and a round pivot, rotating on a smooth floor, for which 
the moment of inertia depends on how the axis tilts. 
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and 

,(8) 



IUA) 



xS aP D^(A) (a, Pel) 



yd aM D { ^{A) 



{a el, (3 ej 

or a e J, Pel) (2.25) 



z5 af} D® (A) + wd^D® (A) (a,[3ej) 
{a = 8 or (3 = 8) 

where 1 = {1, 2, 3} and J = {4, 5, 6, 7}. We denote the SU(3) representation matrix in the 
adjoint (octet) representation as 

D®(A) = ^Tr (A^\ a A\p) (a, (3 = 1, • • • , 8) (2.26) 

and d a 0ry is the usual symmetric tensor. 

The parameters I±, I2, x, y, ~z, and w are to be determined. Note that l\ and I2 are of 
0(1) while x, y, z, and w are of 0(5m). 

The third term in the Lagrangian ()2.21|) comes from the WZW term and gives rise to 
a first-class constraint which selects possible representations. The last term is a potential 
term, 

V(A) = V {l \A) + V {2 \A), (2.27) 
V®{A) = l(l-Dg\A)), (2.28) 

V^\A) = v ( 1 - {d^(A)) 2 ~ (D^(A)) 2 ) , (2.29) 

where 7 is of 0{5m) and v of 0(5m 2 ). Note that V^(A) is of higher order in 5m than 
any other operators. The reason why we include it is that it is the leading order in iV c . 
Equivalently, we assume 

5m 1 

x < m (2 - 30) 

with some relevant mass scale A. 

Collective quantization of the theory is a standard procedure. The only difference comes 
from the fact that the "inertia tensor" I a p depends on the "coordinates" A. The operator 
ordering must be cared about and we adopt the standard one. The kinetic term now involves 
the inverse of the "inertial tensor," and we expand it up to including the terms of order 5m. 
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We obtain the following Hamiltonian, 



H = M d + H + H 1 + H 2 , 



2/i 



2 

a) 1 



(2.31) 
(2.32) 



lf a = xDg>{A)Y,(Faf 



E + E 



da^FaDf^ {A)Ff3 



7=1 



#9 



a&J a,(3£j 7=1 

+ 1 (1 - og'w) , 



where 



2If 2W 
and F a (a = 1, • • • , 8) are the SU(3) generators, 



2/|' 



2/f 



7=1 



1 ' 



^2.33) 
(2.34) 

J2.35) 
(2.36) 



where f a p 7 is the totally anti-symmetric structure constant of SU(3). Note that they act 
on A from the right. 

The WZW term leads to the first-class constraint, giving the auxiliary condition to phys- 
ical states ty(A), 

2 N 

Y R V(A) = -—F 8 ^(A) = -f*(A). (2.37) 

nnnn V3 3 

See m 



52 



53 



54| for more details. In the following, we set N c = 3. 



The novel feature is the existence of the interactions quadratic in F a . Generalizations of 
the Skyrme model have been considered by several authors, but they have never considered 
such complicated interactions as given above. 

It is also important to note that we do not get any interactions linear in F a . The reason 
can be traced back to the fact that the action does not include the terms linear in time 
derivative, except for the WZW term. The absence of such terms comes from the time 
reversal invariance of QCD 7 . 



7 The vacuum angle 8 is assumed to be zero. 
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III. MIXING AMONG REPRESENTATIONS AND THE MASSES OF BARYONS 



A. Symmetric case 

In the absence of the SU(3) flavor symmetry breaking interactions, the eigenstates of 
the collective Hamiltonian furnish the SU(3) representations. The symmetric part H sym = 
M c i + Hq may be written as 

12 



He 



M d + 



1 -)c 2 (SU(2)) + ^-(c 2 (SU(3)) N2 



sym ' 2 \h 

where C 2 (SU{2)) is the spin SU(2) quadratic Casimir operator 

C 2 (SU(2)) = (Faf 



(3.1) 



(3.2) 



with the eigenvalue J( J+l). The operator C 2 (SU (3)) is the flavor SU(3) quadratic Casimir, 



C 2 (SU(3)) = ( F a 



with the eigenvalue 



C 2 ('P, q) = - [p 2 + q 2 +pq + 3(p + q 



(3.3) 



(3.4) 



where (p, q) is the Dynkin index of the representation. Note that we have used the constraint 
F 8 = NJ2V3 in Eq. (JHD- 

The eigenstate is given by the 5*^7(3) representation matrix, 

* M (V; \ J ) ( A ) = V^M^q)(-l) J3 - YR/2 (Y,I,I 3 \D^\A)\Y R ,J,-J 3 y (3.5) 
with Yr — 1, where dim(p, q) is the dimension of representation (p, q) and is given by 

P+ A- (3.6) 



dim(p, q) = (p + l)(q 



1 + 



For the properties of the wave function ()3.5|) . see Appendix 1X1 

By using them, we readily calculate the symmetric mass M-ji of the representation 1Z, 



Mr 



M 



10 



A% = M d + - 



M rI + - 



M r l + - 



5 2 

h + T 2 

3 26 
h + T 2 



3 r i 2 ' 

8 [h + T 2 

3 r i e' 

h + T 2 

"15 14" 
Ji + ~h 

and so on. Because 27 contains the spin-| part and the spin-| part, the former is denoted 
as 27rf, the latter as 21 q . 



M- 



27„ 



M d + - 



(3.7) 
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B. Matrix elements of the symmetry breaking operators 



The SU(3) flavor symmetry breaking interactions mix the representations and an eigen- 
state of the full Hamiltonian H is a linear combination of the (infinitely many) states of 
representations. In this paper, we consider the perturbative expansion up to including 
0(5m 2 ). 

To calculate the perturbative corrections, we need the matrix elements of the symmetry 
breaking operators. Though the calculations are group theoretical, they are considerably 
complicated because of the generators F a . Several useful mathematical tools are summarized 
in Appendix lAl 

Note that the symmetry breaking operators conserve spin SU(2), isospin SU(2) and 
hypercharge U(l) symmetries. So that the matrix elements are classified by (I,Y), the 
magnitude of isospin and the hypercharge, and J, the magnitude of spin. The SU (3) repre- 
sentation and the magnitude of spin are related by the constraint (|2.37[) . 

We introduce the following notation, 



O, = D&\A), O x = Dg{A) (F a ) 2 , (3.8) 

qGX 

8 

Y,d a ^F a D®(A)F p , (3.9) 



E + E 



O z = ^2F a D^(A)F a , O w = E F a Dfrj(A)Fp, (3.10) 

a£j a,f3&J 7=1 

2 



O 



Ul 



and denote them by Oi collectively. 



1. J = \ 

The SU(3) representations with spin | are 8, 10, 27^, and so on. States with the same 
(/, Y) can mix. 

Because all the symmetry breaking operators behave as 8, the octet states can mix with 
10 and 27a in the first order, and also with 35^ and 64^ in the second order. The anti- 
decuplet states can mix with 8, 27^ and 35^ in the first order and also with 64^ and 81 d in 
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the second order. To calculate the masses to second order, we do not need the second order 
mixings, but they are used for the decay amplitudes. 

The matrix elements for the 8 representation are given in Table Note that the matrix 
elements of the operator O z are zero. 



TABLE I: (8\Oi\ 8) 



8 






7 




X 


y 


z 




w 






V2 


N 8 






3 

10 




3 3 
ITS 


20 







VE 

V 

5 




1 
5 


I 

5 


S s 


(1,0) 




1 

10 




3 1 

4 10 


3^ 
20 







3^ 
20 




1 

10 


13 
30 


^8 


(i-1) 




1 

5 




3 1 

4 4 


V3 
5 







20 




1 

10 


1 
2 


As 


(0,0) 




1 

10 




3 1 
4T0 


3^3 
20 







3^ 
20 




1 

To 


3 

To 


Similarly, the matrix elements 


for 10 are g 


iven in Table ITT1 






















TABLE II: 


(lO|0i|lO) 














10 


(I,Y) 




7 




X 


y 


z 




W 




l'l 


V2 


@T0 


(0,+2) 




1 

■1 




3 1 

4 4 


5Vs 

8 


3 
4 




5V3 
8 




1 

7 


3 
14 


Nto 


(s.+i) 




1 

8 




3 1 

48 


5^ 
16 


3 

8 




5V3 
16 




9 

56 


"1 7 
1 t 

56 


^To 


(1,0) 

























1 
7 


8 

21 


^To 


(1,-1) 




1 

8 




3 1 

4 8 


16 


3 
8 




5^ 
16 




5 

56 


25 
56 


In order to calculate the 


masses to second order, we 


also need the matrix elements 


off- 


diagonal in representation. They 


are given in 


Tables ITTT1ITVHV1 andlVri 


















TABLE III: 


(1010,18) 














(I,Y) 


7 






X 


V 


z 




w 




Vl 






(i+i) 


l 


3 1 

4 2^ 


V3 
2VE 


3 

AVE 




V3 

AVE 




1 

AVE 




1 

AVE 


(1,0) 


1 

2VE 


3 1 

4 2V 7 ! 


VL 
2VE 


3 

AVE 




V3 

AVE 









1 

3^ 



Other matrix elements, which are necessary for the calculations of the mixings, are col- 
lected in Appendix [BJ 
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TABLE IV: (27 d |Oi|8) 



(I,Y) 


/ 


x 


■J 


2 


•«! 


Pi 


Do 


(s.+i) 


^6 
10 


3 

4 10 


V2 
5 


4 


^2 
20 


17V3 


13-/3 
140^ 


(1,0) 


1 

5 


3 1 
1 5 


2^ 
15 


1 
2 


30 


1 

70 


9 

70 




10 


3 V6 

4 10 


V2 
5 


4 


V2 
20 


35^ 


_ VI 
7\/2 


(0,0) 


3 
10 


3 3 

4 10 


5 


3 
1 


20 


3 

35 


6 

35 



TABLE V: (27 d \O l \ 10 > 





7 


,r 


y 


Z 


w 


Vl 


V2 


(i+1) 




3 


11 


^3 


23 




V3 




4 8^ 


16\/l0 


2^ 


16Vl0 




56^ 


(1,0) 


1 

4a/5 


3 1 

4 4:VE 


11 

8^ 


1 

VE 


23 


28 


3 

28^ 


(1,-1) 


8^ 


3 

4 8^ 


11 
16V2 


VL 
2V2 


23 
16^ 


56^ 


5%/3 
56^/2 



By using them, we can write down the perturbation theory results for the masses and 



the representation mixings. For example, the nucleon (mainly octet (J, Y) = (|, 
mass is calculated as 



where 



Mn 



M 8 + 77 1 - tt: + -777^ + ^—y - ^—W 



10 



40 



20 



i-i-I 

5 5 



71% - M s 3 M 27j - Mj 



2 Atf 2 V 



27 d 



10-8 



AH k- 



3 x 



3 



3 

To 



V3_ 

4 2v/5 ' 2^5* ' 4^ 4^5 

3 3 v 7 ^ 3 

x y H — 2 H it; 

410 5 4 20 



w 



The masses for other states may be calculated similarly. 



T) state) 



(3.12) 

(3.13) 
(3.14) 



2. J = f 



Spin I states are composed of the representations 10, 27 g , and so on. The matrix elements 
for the 10 representation are given in Tables IVTT1 Note that the matrix elements of O z are 
zero. 
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TABLE VI: (35 d |CM 10> 



(T V\ 


7 


X 


y 


z 


ii J 


V\ 


''2 


(0,+2) 


3 


3 3 

4 4^ 


8V7 


9 

2^ 


8^ 


3 


3 

8^ 


(i+i) 


8^ 


3 \/9 

4 Sv^ 


3V3 

16VT2 


27 
4 V / 14 


3^ 

lev 7 !! 


3 

8VH 


5 

8VTi 


(1,0) 


3 

4^ 


3 3 

4 4^ 


Vs. 


9 

2^/7 


V3_ 

8Vf 


l 

8^ 


11 

24^ 


(1,-1) 


3^ 


3 3^ 




9VE 


VIE 





VE 


8VTi 


4 8^ 


lev 7 !! 


4 V / 14 


16-/L4 


4 V / 14 



TABLE VII: (10|Oj| 10) 



10 




7 


X 


y 


Z 


w 


V\ 


V2 


Aio 


(|,+1) 


l 

8 


15 1 
4 8 


5^ 
16 





VI 
16 


9 
56 


17 

56 


Sio 


(1,0) 

















3 

28 


31 
84 




(i-i) 


1 

8 


15 1 
4 8 


5^ 
16 





16 


5 

56 


25 
56 


n 


(0,-2) 


1 
4 


15 1 
4 4 


5^ 
8 





V3 
8 


3 

28 


15 
28 



The decuplet states can mix with 27 q and 35 in the first order, and also with 35 g , 64 q 
and 81 in the second order. In order to calculate the mass to second order, we need the 
matrix elements of 10 with 27 q and 35 q representations. They are given in Tables IVIIII and 

El 



TABLE VIII: (27 q \Oi\ 10) 





7 


X 


y 


z 


w 


Vl 


V2 


(1,+1) 


VIE 

8V2 


15 VT5 
4 8^ 


5VE 
wV2 


VTE 

8^2 


$VE 

16-^2 


VIE 

28V? 


VIE 

14-v/2 


(1,0) 


1 

4 


15 1 
4 4 


5 

Sv^ 


1 
1 


5 

8^ 


1 

56 


9 

56 


(1,-1) 


V3_ 

8V2 


15 

4 8^ 


5 

16V2 


V3_ 

8V2 


5 

16^ 


V3 
56V2 


5^ 
56^ 



IV. NUMERICAL DETERMINATION OF PARAMETERS 

We can calculate the baryon masses once the Skyrme model parameters are given. In 
the effective theory approach, however, we have to solve in the opposite direction. Namely, 
we need to determine the Skyrme model parameters so as to best fit to the experimental 

1(3 



TABLE IX: (35 1 0*1 10) 





7 


X 


y 


2 


w 


Vl 


V2 


(i.+i) 


5 


15 5 




15 


bV3 


5 


5 




4 8VU 








lev 7 !! 






VE 


15 v 7 ^ 


bVTE 


SVE 


VIE 


VE 


V3E 


(1,0) 




4 4^ 


8Vf 


4V7 


8Vf 




48 






15 3^ 


15-v/IB" 


9VE 







VE 




8VT4 


4 sVTi 


lGv^ 




lev 7 !! 


4v 7 Ii 


(0,-2) 


VE 


15 VE 


5^ 


3VE 


VIE 


VE 


3^ 




iVf 


4 4^ 




AV7 


8Vf 


16V7 


16\/7 



values of the baryon masses. In order to measure how good the fitting is, we introduce the 
evaluation function 

i 1 

where Mj stands for the calculated mass of baryon i, and M^ xp , the corresponding experi- 
mental value. How accurately the experimental values should be considered is measured by 
(7j. Because we neglect the isospin violation effect completely, we use the average among the 
members of an isospin multiplet for the mass and the range of variation within the isospin 
multiplet for the a*. This is why our estimate of cij for isospin singlets is severe, while <tn is 
considerably large though the masses of proton and neutron are very accurately determined. 
At any rate, these numbers should not be taken too seriously. 

The sum is taken over the octet and decuplet baryons, as well as + (154O) and 0(1860). 
Note that we have nine parameters to be determined. We need at least one more state than 
the low-lying octet and the decuplet. In this sense, our effective theory cannot predict the 
@ + mass. In our calculations, we use the values of M^ xp and <Tj given in Table 1X1 

TABLE X: Experimental values of baryon masses and their deviations used in our calculation. In 



the last low, we 


also give 


the baryon 


masses 


calculated by using 


the best fit set of parameters 


621). 


(MeV) 


N 


£ 




A 


A 


£* 




n 


e 




M exp 


939 


1193 


1318 


1116 


1232 


1385 


1533 


1672 


1539 


1862 


<*i 


0.6 


4.0 


3.2 


0.01 


2.0 


2.2 


1.6 


0.3 


1.6 


2.0 


Mi 


940 


1180 


1332 


1116 


1228 


1389 


1537 


1672 


1539 


1862 



The problem is a multidimensional minimization of the function \ 2 of nine variables. In 
general such a problem is very difficult, but in our case, thanks to the fact that the function 
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is a polynomial of the variable in the perturbation theory, a stable numerical solution can be 
obtained. Our method is basically Powell's one, but we tried several minimization algorithms 
with the equivalent results. The best fit set of parameters is the bottom point of a very 
shallow (and narrow) "valley" of the function, and x 2 does not change very much even if we 
change the values of parameters in a certain way. 
The best fit set of parameters is 

M d = 389 MeV, if 1 = 174 MeV, Jj 1 = 585 MeV, 7 = 832 MeV, 
x = 27.8 MeV, y = -104 MeV, z = -306 MeV, w = 111 MeV, 
v = -148 MeV, (4.2) 

and leads to the mass spectrum given in Table IXl with \ 2 — 4-4 x 10 1 . 

Note that best fit set of parameters is quite reasonable, though we do not impose any 
constraint that the higher order (in 5m) parameters should be small. The parameter 7 is 
unexpectedly large (even though it is of leading order in N c ), but considerably smaller than 
the value (7 = 1573 MeV) for the case z seems 

also too large and we do not know the reason. 

Once we determine the best fit set of parameters, we can predicts the masses of (mainly) 
anti-decuplet members, 

M N < = 1711 MeV, M E / = 1819 MeV. (4.3) 

Compare with the xQSM prediction Q], 

M w = 1646 MeV, = 1754 MeV. (4.4) 

It is tempting to identify N' with N(1710) (* * *) state. On the other hand, for E' there 
are two candidates, £(1779) (*) and £(1880) (**). In any case, it is too early to identify 
them. 

The mixing coefficients of the eigenstates are also obtained. For the (mainly) octet states, 
the coefficients are given in Table IXTl 

Similarly, for (mainly) decuplet and (mainly) anti-decuplet states, they are given in Ta- 
bles EH] and |Xffl] The numbers in parentheses are the second order contributions calculated 
by using the parameters given in (J4.2)) , which are used in Sec IV Bl 



18 



TABLE XI: Mixing coefficients for the (mainly) octet states in the first order perturbation theory. 
The numbers in parentheses are those in the second order. 





N 


E 




A 


8 


1 (-0.056) 


1 (-0.051) 


1 (-0.014) 


1 (-0.021) 


10 


0.288 (0.036) 


0.288 (0.086) 


0(0) 


0(0) 


27 d 


0.169 (0.087) 


0.138 (0.087) 


0.169 (0.036) 


0.207 (0.075) 


35d 


(0.085) 


(0.080) 


0(0) 


0(0) 




(0.038) 


(0.028) 


(0.038) 


(0.051) 


TABLE XII: Mixing coefficients 


for the (mainly) decuplet states in the first order perturbation 


theory. The numbers in parentheses are those in the second order. 




Ui 


A 


E* 




n 


10 


1 (-0.150) 


1 (-0.082) 


1 (-0.033) 


1 (-0.003) 


27 9 


0.543 (0.122) 


0.397 (0.031) 


0.243 (-0.017) 


0(0) 


35 


0.065 (0.076) 


0.082 (0.072) 


0.087 (0.051) 


0.082 (0.024) 


35 g 


(0.138) 


(0.049) 


0(0) 


0(0) 


64, 


(0.053) 


(0.050) 


(0.033) 


0(0) 


81 


(0.009) 


(0.013) 


(0.015) 


(0.013) 


The 


mixings are rather lar£ 


;e. One may think that the perturbation 


theory does not 



work. For the (mainly) octet states, the second order contributions are much smaller than 
the first order ones, while for the (mainly) decuplet and the (mainly) anti-decuplet states, 
the mixing with 27^ and 35^ are large and the magnitude of the second order contributions 
are comparable with the first order ones. 

V. DECAY WIDTHS 

In this section, we calculate the decay widths of various channels based on the calculations 
done in the previous sections. Since our treatment of the baryons is a quantum-mechanical 
one, the full-fledged field theoretical calculation is impossible. What we actually do is a per- 
turbative evaluation of the decay operators in the collective coordinate quantum mechanics. 
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TABLE XIII: Mixing coefficients for the (mainly) anti-decuplet states in the first order perturba- 
tion theory. The numbers in parentheses are those in the second order. 





e 


N' 


£' 





8 


0(0) 


-0.288 (-0.089) 


-0.288 (-0.157) 


0(0) 


TO 


1 (-0.061) 


1 (-0.160) 


1 (-0.235) 


1 (-0.285) 


27 d 


0(0) 


0.314 (0.304) 


0.513 (0.219) 


0.703 (-0.079) 


35d 


0.351 (0.130) 


0.372 (0.193) 


0.352 (0.233) 


0.277 (0.225) 


64 d 


0(0) 


(0.136) 


(0.203) 


(0.215) 


81 d 


(0.116) 


(0.128) 


(0.116) 


(0.082) 



A. Formula for the decay width 



Since there seem to be confusions |55l. l56l l57| concerning the factors in the decay widths, 
we reconsider the derivation of the formula. See ^| for the discussions of the calculations 
of the decay widths in the Skyrme model. 

Decay of a baryon to another baryon with a pseudoscalar meson may be described by the 
interactions of the type 

£>decay = -igd^ , (5.1) 

where is the baryon axial-vector current and Q is a pseudoscalar meson field. The 
coupling g has the dimension (mass)~ l and usually related to the pion decay constant F n , 
g ~ F~ l . In the nonrelativistic limit, the time component may be dropped, and it is useful 
for us to work in the Hamiltonian formulation, 

H decay (t) = j d 3 xH decay =ig J d 3 xd k r(x)jf(x). (5.2) 

In the leading order, the amplitude of the decay B — > B'cf) may be given by 

A = J dt{<t> a (p)B'(P')\H decay \B{P)) 

= ig J d 4 x(cf> a (p) \dk<P{x)\ 0) (B'(P') \j/ k (x)\ B{P)) . (5.3) 

With the relativistic normalization of the state, the first matrix element may be written as 

(<j) a (p) \d h <jp(x) | 0) = ip k e ipx 5 al3 . (5.4) 
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In our treatment of baryons, the state \B(P)) has the wave function 



^2E B {P)m B {A)e~ iEBi - p)t+%px , (5.5) 

where E B (P) = \/ P 2 + M B and X is the position of the baryon. The state satisfies the 
relativistic normalization, 

(B'(P')\B(P)) = 2E B (P)(2n) 3 5 3 (P' - P), (5.6) 

where the inner product is defined as the integration over X and A. The axial- vector current 
may be obtained from the action by using Noether's method. After replacing U(x) with 
A(t)U c (x — X)A^(t), we obtain the collective coordinate quantum mechanical operator by 
representing the "angular velocity" u a in terms of the generator F a ("angular momentum") 
according to the usual rule. J£ k depends on F a but not on t anymore. Note that it depends 
on X through the combination x — X. Now the second matrix element may be written as 

{B\P')\j^\x-X)\B(P)) 
= ^2E B ,2E B e t{E B'- E ^ t j dA** B ,(A) d 3 Xe^ p '~ p > x J b pk {x - X) j V B (A).(5.7) 

By making a shift x — »• x + X and integrating over X and t, we get 8 

A = (2vr) 4 5 4 (P / + p-P) y / 2E B ,2E B J dAV* B ,(A) (gp k j d 3 xe- ip - x J 5 ak (x)^j V B (A) 
= (27r) 4 5 4 (P' +p- P)M. (5.8) 

In the nonrelativistic approximation, we put e~ tp x ~ 1. The integral of the axial- vector 
current, 

jO% ecay = gp k J d 3 xJ 5 ak (x), (5.9) 

depends on A and F a , and has the right transformation property under flavor SU(3) trans- 
formations. Here we introduce a mass scale A, which we take A = lGeV. It is known [so| 
that the leading order result is given by 

jOl cay ^gCD^ k (A)p k , (5.10) 



M. corresponds to the invariant amplitude in the relativistic field theory, though, in our formulation, 
relativistic property has already been lost. The spin of the wave function (|5.5(l docs not transform 
properly under boosts. The decoupling of spin reflects the fact that baryons are now considered to be 
(almost) static, i.e., \P\ < M B and \P'\ < M B >- 
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where C is a dimensionless constant. In the "traditional" approach, this is a functional of 
the profile function F(r), and obtained by explicitly integrating J£ k (expressed in terms of 
U c (x), A, and F a ) over x. Note that it has no dependence on Mb nor on M^. 

In the effective theory approach, on the other hand, the coefficient of the decay operator 
may be determined by fitting the widths to the experimental values. We follow this way. 

The leading operators are well-known, 



decay 



G D®(A)-G 1 d kM D%(A)F,-^Dtl(A)F t 

v 3 



V 



(5.11) 



where the index k runs 1,2,3. The constants G a (a = 0, 1,2) are dimensionless, thanks to 
the explicit mass scale A. 

Once we calculate the amplitude /A, we are ready to obtain the decay width, 

1 3 



IPI 



8ttM b A2 



1 >M? 



(5.12) 



where \p\ stands for the magnitude of the meson momentum in the rest frame of the initial 
baryon B, 

-L + M b>) 2 - m 2 ] [(M B - M B ,y - m% (5.13) 

ZlViB 



where m is the mass of the meson. \M.\ is defined as 



\M\ = ; 

\P\ 



(5.14) 



sptn,isospin 



The symbol 



spm,isospm 



denotes the average of the spin and the isospin for the initial state 



baryon as well as the sum for the final state baryon. By extracting the normalization factor, 
4M B E B >(\p\), 



\M\ =AM B E B >{\p\) 



M 



we may rewrite it as 



B^B' 



\pf E B j\p\) 
2tt A 2 M b 



(5.15) 



(5.16) 



This is our formula for the decay width. 
The widely used formula, 



B^B' 



A-kMbMb 



(5.17) 
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which corresponds to Eq. ()5.16|) seems to be based on the interaction of Yukawa type 



C = g<pB'B^ B 'h5^ a ^B(p a - (5.18) 

It is this coupling constant g^s'B that depends on the initial and final baryons. Actually, 
Goldberger-Treiman relation relates g^B'B wit ha in Eq. (j5.1j) . From this point of view the 
inverse mass factors in the coefficients in Ref. [2( may be understood. Most of the authors 
do not seem to care about the normalization (relativistic or nonrelativistic) of the states. 

We have a preference to Eq. (|5.1j) over Eq. (j5.18|) because the derivative coupling is a 
general sequence of the emission or absorption of a Nambu-Goldstone boson at low-energies. 
On the other hand, the universality (i.e., independence of the initial and final baryons) of 
the coupling g is generally less transparent. It is however very naturally understood in 
the Skyrme model, in which, as we see above, the axial-vector current comes from a single 
expression for all the baryons. 

It is interesting to note that, although the reasoning seems very different from ours, the 
decay width formula with the ratio of baryon masses in Ref. Q], 

|P|S Mb ''M\ 2 (5-19) 



" 2ir(M B + M B ,f M B 

looks similar to ours (|5.1fij) . if the factor (Mb + Mb>) is identified with our common mass 
scale A. 



B. The best fit values of the couplings and the predictions 



In this subsection, we calculate the important factor 



M 



and then the decay widths. 



The calculation goes as follows. What we need to calculate is the matrix element, 



(V B > \O a decay \ tf B ) = J dA^^O^! 



{A). 



(5.20) 



The baryon wave function ^b{A) is a linear combination of the states in various represen- 
tations, 

*b{A) = (5.21) 

i 

where V%{A) is the eigenstate of H with F = (Y, I 3 , 1) and S = (Y R = +1, J 3 , J) in 
representation TZi and the coefficients cf is those we obtained in the previous section. So we 
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first calculate (^fpf S i \02 ecay \ ^"fs/- Furthermore, the spin and flavor structure is completely 
determined by the Clebsch-Gordan (CG) coefficients. For example, consider the matrix 
elements of the G decay operator, 



dA(^ s ,(A)) D%(A) P ^ S (A) 



(8), 



dimRi 
dim^ 



P*(S')P(S)J2( I 



F F' I \ k S S' 

r 



V 



(5.22) 



Because 



8 Hi TLj 
k S S' 



V 



1 ■ 2 
P — ip 



V2 



(J',-J 3 \1,1;J,-J 3 ) 



p 1 + ip 2 



(J',-4\1,-1;J,-J 3 ) 



V2 

-p 3 (J',-4|l,0; J,-J 3 ) 





( 8 


Hi 










V 10 


J 1 


J,1 Jr 



(5.23) 



where (J, J 3 \j, j 3 ; j', j' 3 ) is the usual SU(2) CG coefficient and the last factor is the SU(3) 
isoscalar factor, the spin factor [■ • • ] can be extracted. In this way, the matrix element may 
be written as 



dA [^ s ,(A)j O a decay mf s {A) = (isospin CG) x [spin CG with p k ] x AQ^, (5.24) 

where •M a \p % i p\ contains all the other factors, such as the phases, the isoscalar factors, and 
so on. The first two factors are irrelevant when we calculate the average and sum of the 
spins and isospins and just give | |p| 2 . Only M^ 3 ^) ^ matters. Actually, after squaring the 
amplitude, averaging the spin and the isospin of the initial baryon, and summing over the 



M 



as 



spin and isospin of the final baryon, we have 

£( C fr^> )C f 



M 



2 _ 1 

" 3 



(5.25) 



In the following, we give the factors A^ Q | F , 1 F \ for various decays in matrix forms. Before 
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presenting the factors, we introduce the following combinations of couplings 



Gs,10 = 


= G + 


2 Gl ' 


Gs,27 = 


- Go - 


1 /""t 

- 5 Gi, 


^10,27 


— Go 


+ G X . 


^10,35 


= Gq - 


- G\, 


^27,10 


= Gq 


-2d, 


Lt 27,27 


= G\ 




^27,27 


— Gq - 


' 


^27,35 


— Go 




^27,35 


— Go 


-2G X 


^27,64 


— Gq - 


■ 


^35,10 


— Go 


— Gi, 


^35,27 


— Go 


-2G X 


r (l) 
35,35 


= G -f 


■HGi, 


r (2) 
35,35 


= Gq 


+ 5^1, 


^35,64 


— Go 


+ 2 Gl 


^64,27 


= Gq - 


■ 


^64,35 


= Gq 


5 /^r 
- g^l, 


^64,35 


— Go 


+ 1 G 1 


^64,64 


— G\, 




^(2) 
^64,64 


— Go 


- 2^1) 


^64,81 


— Go 


+ f^l 



^10,8 1 


- Go - 




-^10,27 


= Go 


-2G l5 


-^27,8 1 


- Go - 


- |Gi, 


-^27,10 


= Gq 


+ G U 


r 27,27 


= Gi, 




P (2) 
r 27,27 


= Go 


- 5G1, 


-^27,35 


= Gq 


-2d, 


-^27,64 


= Go 


- |Gi, 


-^35,27 


= Go 


+ G 1; 


-^35,64 


= Gq 


- 5^1, 


-^35,10 


= Go 


— Gi, 


-^35,27 


= Go 


-2G?i, 


F w 

35,35 


— Go 


+ HGi, 


F (2) 
35,35 


= Go 


+ |Gi, 


-^35,64 


= Go 


+ §Gi, 


-^35,81 


— Gq 


— G\, 


-^64,27 


— Go 


2 ' 


-^64,35 


= Go 


+ §Gi, 


r 64,64 


— G\. 




F (2) 
r 64,64 


= Gq 


-^64,81 


= Gq 


5 /^r 


-^81, 64 


= Go 


+ |Gi, 















(5.27) 



^8,8 " 


= G + 


\Gi - 


|g 2 , 


#(2) _ 
-"8,8 - 


= G 4 


- 2G1 + 


|g 2 , 


-^8,10 = 


= Go - 


-^8,27 1 


- Go - 


-2Gx + 


|G 2 , 


-^10,8, 


= Go 


-Gi- 


|g 2 , 


-^10,10 


= Go 


-^10,27 


= Go 


+ ^1 


n G 2, 


-^10,35 


= Go 


+ |Gi - 


- |g 2 , 


-^27,8 = 


= Go - 


-^27,10 


= Go 


+ n G i 




-"27,27 


= Go 


, 13,-y 
+ "g-Gi 


+ |G2, 


-"27,27 


— Go 


-"27,35 


= Go 


17 ri 
- Yo G l 


3 p 
~ T0G2, 


-^27,64 


= Go 


-2Gi- 


f- f G2, 


^35,10 


= Go 


-^35,27 


= Go 


- Yo^i 


- ^G 2 , 


rr(l) 

35,35 


= Go 


11 ,0 

- T Gi 


+ |G2, 


rr(2) 

35,35 


= Go 


-"35,64 


= Go 


+ |Gi - 


f |g 2 , 


-^35,81 


= Go 


+ |Gi - 


- |g 2 , 




= Go 


-^64,35 


= Go 


+ |Gi - 


f |g 2 , 


#(1) 
-"64,64 


= Go 


+ f G l " 


f- ^G2, 


-"64,64 


= Go 


-^64,81 


= Go 


16 r< 
— 7-^1 


3 r 















Gi — 4G2, 



- |Gi + |G2, 
2Gi + |G2, 

_ 29/^< _ 13/o 
38 1 38 2) 



\G, - |G 2 , (5.28) 

1 

2 

G2, 



19 r< _i_ i/^ 

22^1 + 2^2 
2Gi+| 

19^ 9_^ 

22 1 22 2 



which are useful in representing M^^] F y Our naming conventions: G for the decay of the 
(mainly) decuplet to the (mainly) octet, F for the decay of the (mainly) anti-decuplet to 
the (mainly) decuplet, and H for the decay of the (mainly) anti-decuplet to the (mainly) 
octet. The subscript implies the components. The superscript (1) or (2) distinguishes the 
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outer degeneracy. For example, Gs,io stands for the coupling which appears in the matrix 
elements between 10 and 8, and is needed for the calculation of the decay of a (mainly) 
decuplet baryon to a (mainly) octet. 



The second order results, which 
59|. Our results are very extended and 



The first order results have been given in Refs. 
overlap with some of ours have been given in Ref. 
lengthy. But they are actually very important part of the present paper, and useful for 
xQSM calculations too. 

For the decays of (mainly) decuplet baryons, we have 



A -> N 



- 7T 



/ 3 r 



V 










^10,27 

3^(2) 
7^27,27 

i r 

2v / l05 Cj35 > 27 






5 n 

2 721 G27 > 35 



3 p 

^^10,35 

Vh ri _ 
2V21 27 > 35 
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Go 



J_nW_ , 25/^(2) zVlfi 
68^35.35 34^35.35 4^ LT 35; 



rv7 G64 > 35 



35^64,35 



4^(2) 
7^64,64 



\ 





1^64,81 / 



(5.29) 



£* -> A + 71 : 



( 



v4o Gs > 10 

V2 ft 



\ 







V2 r 











^7g^27,27 V? G 27,35 

1 r~i 3\/5 r~i 3 r~i 

5g < -J r 64,27 ir-m^PA?^ — Lr, 



27,35 
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G 



27,64 



^^64,35 ^^64,35 



42V 7 ! 
2^5^(2) 
7 ^64,64 








64,81 J 



7T : 



7^8,10 



275 G27 ' 10 










2 ^10,27 
4 U 27,27 4y /j 



<JQ7 97 -r7S<-J27,35 



4^/7^35,27 




2 V / 14 



G 



64,35 





2^/7^10,35 



_V15^ _ 
4V7 27 > 35 









1 _i__5_r»( 2 ) VT5^_ 

' Sv^ 35,35 4^ 35,35 2^7 35 > 64 



2^ 64 > 35 



75 










3^ 



^5^64,64 2v ^q 



^64,81 / 



(5.31) 
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-> H + tt : 

-1^^8,10 



10 

_j_n 57 ^(2) 

ITTs^ 27 - 10 WTo^ 27 - 27 " 







15 nt 1 ) 25 n 
Vl3 G 27,27 47 42 G27 ' 35 

9 





5^5 q 



28^ 



27,64 









28^ 



64,27 



4 V / 14 



r< -3 r^ 1 ) i ii r<(2) 7^ ^ 

fi4 > 35 2V5 ^64,64+ 7^ ^64,64 



64,81 



/ 



With the mixing coefficients, one can easily calculate 
A — > N + tt, it is given to first order by 



M 



„N „N N 
c 8 ^0 C 27 d 



/ 3 r> V2, 




^10,27 



M 








(5.32) 

For example, for the decay 



\ 



1 n 3/^(2) 5 r< 

\730 G27 ' 10 7 G 27,27 2 73T G27 ' 35 



10 



"27„ 



V W 



(5.33) 



with Cj and c~ being given in Tables IXTI and IXTTl Note that the size of the matrix depends 
on the quantum numbers (I, Y) of the initial and final baryons, and corresponds to our 
mixing coefficients given in Tables IXllXlll an dXilll 

The factors for the (mainly) anti-decuplet baryons are obtained similarly. First, for @ + , 



+ 



-> N + K : 

3 TT_ 

'4-" 10, 10 
7V3 



H 



4VW "27,10 
1 n 



10 



V 



4^14^35,10 




3 TT_ 

~4V7 ^10,35 

5^ TT 

' 4v ^4^27,35 
(1) 
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17%/2 35,35 68^ 35,35 
V70 ^64,35 







yj.i7 

^T-" 35,81 



7 TT 

'V30 ^64,81 



(5.34) 



There are several interesting decay channels for the (mainly) anti-decuplet excited nucleon 
N', for the (mainly) anti-decuplet and for the (mainly) anti-decuplet <f> decays, 

Let us assemble all of these ingredients. In perturbation theory, it is important to keep 
the order of the expansion. In order to obtain the widths to second order in 5m, we need to 
drop the higher order terms consistently. 

Because the mixings among the representations are large, however, we often get negative 
decay widths if we drop higher order terms in the squares of the amplitudes. Of course such 
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N' -> N + vr : 

/ _27 rr(l) _ 3rr(2) 
20""8,8 4- fl 8,8 
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2^^10,8 

"sTS^ 27 ' 8 
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8 10,10 
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49 TT_ 

8^30^10,27 
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9 H. 
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25^ TT 

W2l 27,35 
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£' -> N + K : 
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results are unacceptable. Considering that the perturbative contributions are large, we also 
show the results in which the squares are not expanded. (In Appendix [0 we show that 
negative decay widths can appear even the mixings are smaller.) 

We perform the second order calculation as well as the first order one. In the following, 
the case (a) is the first order result with the square of the amplitude being expanded, the 
case (a') with it being not expanded. The cases (b) and (b') are corresponding second order 
results. The second order results are drastically different from the first order ones. But, 
before discussing the results, let us explain the procedure. 

First of all, we best fit the couplings, Go and G%, to the experimental decay widths of 
the (mainly) decuplet baryons. The procedure is similar to that in the previous section. We 
minimize 

/■p r exp\2 

* 2 = E ( J' } ■ < 5 - 5 °> 



where Tj stands for the calculated value of channel i, T^ xp is its experimental value, and Oi 
represents experimental uncertainty. The values we use are given in Table IXIV1 We give 
the best fit set of couplings Go and G\ in Table IX VI We then determine Gi by using the 
F/D ratio, 

F/D= 6 -^. (5.51) 

-"88 

The experimental value of F/D is 0.56±0.02, though many authors prefer to use F/D = 0.59. 
The values of G2 for various cases are also given in Table IXV1 

TABLE XIV: Experimental values of the widths of the (mainly) decuplet baryons and their uncer- 
tainty. 



(MeV) 


A -> N + 7T 


S* -» A + vr 


S* -» S + 7T 






120 


32.6 


4.45 


9.50 


a 


5 


0.74 


0.74 


1.8 



The best fit values of the (mainly) decuplet decay widths are summarized in Table IXVI| 
where we also give the phase space factor K, 



r = k 



M 



lp| 3 E 



Ks %»w B - (5 - 52) 
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TABLE XV: Coefficients of the decay operators. Go and G\ are obtained by fitting the calculated 
decay widths of the (mainly) decuplet baryons to the experimental values, while G2 is fixed by the 
F/D ratio. The results (a) and (a') are those in the first order with the squares of the amplitudes 
being expanded, and not expanded, respectively. The results (6) and (b') are those in the second 
order. The x 2 values for these cases are also given. 





(a) 


(«') 


(b) 


(&') 


Go 


4.74 


5.29 


5.64 


5.77 




13.5 


8.81 


15.4 


9.93 


x 2 


38.8 


4.73 


13.1 


6.43 


G 2 


0.07 


0.08 


0.08 


0.09 



We use the experimental values for the baryon and the meson masses, if they are known, 
in all decay width calculations. For (mainly) anti-decuplet N' and we use our predicted 
values (g3J). 



TABLE XVI: Decay widths for the (mainly) decuplet baryons. Calculations are done with the 
coupling constants given in Table IXVl The kinematical factor K is explained in Eq. 1)5.52(1 . 



(MeV) 


K 


r (o ) 




r (6 ) 


I>) 


A -> Nvr 


1.47 


92.1 


114 


110 


112 


£* -> Avr 


1.18 


33.9 


32.8 


32.8 


32.9 


X* -> Svr 


0.26 


5.08 


5.50 


5.73 


5.31 




0.49 


13.0 


11.4 


13.9 


12.2 



Once we determine the best fit values of G a , the decay widths of (mainly) anti-decuplet 
baryons can be calculated. The results are given in Table IXVIII 

By looking at the Table IXVlIl one can easily see that the widths change rather randomly 
by going to the second order from the first order. It gets wider in a channel, while narrower 
in another. The behavior also depends on whether we expand the amplitude or not. 

We have also done the calculations with the calculated baryon masses even for N, A, etc., 
and have found that the results change drastically. 

Our result for the width of B + is of order 100 MeV, which clearly contradicts with 
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TABLE XVII: Predictions for the decay widths for the (mainly) anti-decuplet baryons. 



(MeV) 


A 
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i (*0 


1(6) 


L (b') 


e+ . 
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21 / 


y A 1 
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141 ( 


1U6 
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£i VI 


A 1 
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U.U 7 


Q A ^ 
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U.61 


N' - 


-> A7T 


( .36 


U 


vo a 
( 2.6 
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N' - 
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Z.U4 
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zo.l 


Q Q7 
O.O f 


5o.5 


Q A d 

o4.0 
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n 00 
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2.5U 
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r in 

5.1U 
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1 r 1 

15.1 


1 ^7/1 1 

-L(4l 


zz.l 


1 ^ /i 
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85. / 


E' - 


-> Evr 


14.0 


A AO 

44o 


86.6 


to ( 


5U. / 


E' - 


-> E77 


1.62 


/I 1 1 

-4.11 


n £M 
U.61 


A O A 

-13.4 


n 1 9 
U.13 


E' - 


-> A7T 


18. / 


A Id 

4/d 


114 


1 21 


1 n A 
ly4 


h - 


1 — 'TV 

-> aK 


n HQ 
U.Uo 


U.1U 


U.U4 


n 1 1 
U.ll 


n n/i 
U.U4 


E' - 


+ £*7T 


5 94 


n 


6 71 


2.25 


81 1 


£' - 


■+ AK 


1.93 





41.2 


25.2 


116 


4>- 


* SK 


5.12 


-137 


14.9 


-181 


44.5 


<\>- 


■» H7T 


10.8 


416 


62.9 


294 


58.6 


<\>- 


> a 7r 


2.68 





8.63 


47.2 


12.5 



reported experimental results. 

Still, one may get some insights from our calculations. First of all, we see the smaller the 
coupling G\ is, the narrower the widths are. Second, our results are almost insensitive to 
the value of G 2 . We did the similar calculations with F/D = 0.59, but the results are only 
slightly changed and qualitatively the same despite the fact that the value of G 2 changes 
considerably. 

VI. SUMMARY AND DISCUSSIONS 

In this paper, we reconsider the Skyrme model from an effective theory point of view. 
In this approach, the Skyrme model parameters, which appear in the collective coordinate 
quantized Hamiltonian, are determined by fitting the calculated baryon masses to the ex- 
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perimental values. Once the Skyrme model parameters have been fixed, various physical 
quantities can be calculated. In particular, we make a prediction for the masses of N' and 

The idea behind this approach is that the %PT provides a framework that represents QCD 
at low energies in which the soliton picture of baryons, which is a consequence of the large- N c 
limit, emerges. We start with the action up to 0(p A ) and keep only the terms which are of 
leading order in N c . After quantization, we calculate everything as a systematic expansion 
in powers of 8m. Note however that we keep in mind that there are infinitely many terms 
which contribute to the Skyrme model parameters. Thus the number of parameters appears 
to increase from the starting action to the Hamiltonian. From the effective theory point 
of view, it is not the number of independent parameters of the "model" but the symmetry 
and the power counting that matter. The "derivation" given in Appendix is, however, 
convenient to generate relevant operators which respect them. 

The basic idea that the higher order contributions improve the Skyrme model picture 
seem to be justified in Appendix O by comparing the conventional Skyrme model with the 
one with L5 and Lg terms in the "traditional" approach. 

We have performed the complete second-order calculations for the masses, and determined 
the Skyrme model parameters. We find that, although the octet behaves good, the decuplet 
and the anti-decuplet have large mixings and the perturbative treatment may be questioned. 

We also re-examine the decay calculation in the Skyrme models, by deriving the formula 
for the decay width from the derivative coupling interaction to Nambu-Goldstone bosons. 
A careful derivation reveals how the decay widths depend on the initial and final baryons. 

We calculate the widths of several interesting decays by using the formula. In particular, 
our calculation predicts a wide decay width for + , in contradiction to the experiments. If 
@ + has really a very narrow width, as reported, our theory fails to reproduce it. A possible 
explanation of this failure is that our perturbative treatment is poor for the (mainly) decuplet 
and the (mainly) anti-decuplet states. Because the decay parameters are determined by using 
the (mainly) decuplets, this could influence very much. Another explanation comes from the 
very subtle nature of the decay width calculations. The results heavily on the kinematics, 
i.e., the masses of the baryons and the factors in the formula. A few percent change of 
the mass can often cause a hundred percent (even more) change of the decay width. The 
theoretical ambiguity is extremely large. The results in Appendix O seem to support this 
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explanation. 

Is our fitting procedure appropriate? We vary all of the parameters as free parameters 
and treat then equally. But there must be a natural hierarchy in them: leading order 
parameters must be fitted to the bulk structure, and subleading parameters should account 
for fine structures. We tried to find such a systematic procedure, but so far, the presented 
method is the most satisfactory. 

What should we do to improve the results? It is difficult to go to the next order in 
perturbation theory, because in the next order we need to include more operators, thus 
more parameters to be fitted. A diagonalization, rather than perturbative expansion, may 
be an option, though it somehow goes beyond the controlled effective theory framework. 

Does + really exist? Is it narrow? Why so? Is the narrowness a general feature of 
the Skyrme model? We do not have definite answers yet. But our results suggest that if it 
really exists and is really narrow, it seems very peculiar even from the Skyrme model point 
of view. As shown in Appendix El it is not just because of the difference of the symmetry 
breaking interactions. 



Praszalowicz |24j discussed how G + becomes narrow in the large- N c limit, and showed that 
the narrowness comes from the interplay between the cancellation in Gjq = H 8 jq (which 
becomes exact in the nonrelativistic limit) and the phase space volume dependence. The 
important factor of his argument is of course the cancellation in Gjq, but it comes from the 
xQSM calculations. In our effective theory treatment, on the other hand, the couplings G a 
are parameters to be fitted. The only possible way to understand such "cancellation" in the 
effective theory context seems symmetry. We do not know if it exists, nor what it is. 

APPENDIX A: MATHEMATICAL TOOLS 

In this section, we summarize some basic mathematical formulae for the calculations of 
the matrix elements. 

1. Properties of the wave function and basic formulae for matrix elements 

Let us first introduce the notations. A baryon wave function \1/ has the flavor index 
F = (Y, I3, 7) and the "spin" index S = (Yr, J3, J). The eigenstate wave function ()3.5j) of 
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Hq may be written as 

V n (g) FS = VdhtflP(S) (p^jsig))* , 9 e SU(3), (Al) 

where P(S) = (— 1) J3_Yr ' 2 , S = (Yr, —J3, J), and D n (g) is the representation matrix of g 
for representation 7Z. Note that physical states must satisfy the constraint ()2.37[) . Yr — 1, 
but this is irrelevant to most of the results in this section. The wave function is normalized 
in the sense, 

(Vfs^f's') = <W<W<W (A2) 
Flavor transformation / G SU (3) acts from the left, g 1— > fg. The corresponding unitary 
operator Ufi avor (f) acts as 

U flavor (m n (g) FS = ^ n {r l g) FS = y n {g)F>sDl F {f). (A3) 

Here and hereafter, the summation over repeated indices is understood. On the other hand, 
"spin" transformation 9 s G SU(3) acts from the right, g > gs" 1 . The corresponding unitary 
operator U spin (s) acts as 

U spm (s)^ n (g)Fs = V n (gs) FS = Vdh^fZP(S) (D* 8 ,(g)y {Df, § (s))* 

= Vd^KP(S') (D^g))* Df,* s *(s), (A4) 

where 7Z stands for the conjugate representation to 7Z and S* stands for (—Yr, J 3 , J). We 
have used the phase convention of Ref. j^. When s is restricted to the "upper-left" SU(2) 
subgroup, it reduces to the usual spin transformation law. 

The infinitesimal transformation (Lie derivative) of the "spin" transformation defines the 
operator F a introduced in Eq. (|2.3b|) . 

F a ^(g) FS = v / dh^JZP(S) (D%(g))* (-T?) St§ , , (A5) 
where is the SU (3) generator in the representation 7Z. In particular, because 

we have 

F^ n (g) FS = -^Y R * n (g) FS . (A7) 



Only the SU (2) subgroup corresponds to the usual spatial rotation. 
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The basic calculational tool for various matrix elements is the orthogonality of irreducible 
representations, 

J dg {D*{g)y DSig) = -^5 n , w 8 ik 5 jU (A8) 

where dg is a normalized Haar measure. For a compact group such as SU (3) it is left- and 
right-invariant. 
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, |6J. It 



Another important tool is the ST/ (3) Clebsch-Gordan (CG) coefficients jool. 
enables us to calculate the following integral, 

/ Mm^'^MM = a^g (: : :)] (: : :) r . m 

where the subscript r counts the multiplicity m of the representation 1Z in the direct product 
representation 1Z\ ® IZ2, or, equivalently, the multiplicity of IZi in 1Z <8> 7^2- Note that in 
the following we do not always work with the "physical basis" which diagonalizes the (right- 
)hypercharge and (iso-)spin, the CG coefficients are not necessarily real. 

All of the operators whose matrix elements we need to evaluate involve the octet (adjoint) 
representation. Thus our first formula is 



* FS 



D (8) 



<s) = Jdg {*%M) m D%)*?: si (g) 



Vd^VdhtiR>P*(S)P(S') ( / dg {D%(g)Y D%{g)D%{g) 



rp 2. ^ / T* \ / T 

A decay operator contains (at most) an with thus our second formula is 



v / d^v / dim^P*(5)P(5') (| dg (D^))* D^D*^)) (-7f' ).,-„ 



dimft ' ^ \ a F' F I \ (3 S" S I \ / S',S 

r=l x 7 ' 

Most of the symmetry breaking operators in ifi contains two F a 's with P/ 8 -*. The third 
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formula is useful in evaluating the matrix elements of them, 



n 

FS 



F S D a/3 F V 



^F'S' 



VdimKVdimR'P* (S)P(S') 



S',S[ 



dim7?. 



/8 5; Si 



T, 



R' 



S',SJ 



The sum of the last three factor may be rewritten as 



(If): 



8 TV TZ 



T, 



s',s; 



(A12) 



(A13) 



2. Simplification in the diagonal case 



Certain simplification occurs in the case TZ = TZ'. In this case the CG coefficients may 
be considered as dim7?. x dim7?. matrices, 




(A14) 



In the following, we derive several useful formulae by giving the explicit forms of the matrices 
Ma ■ The analysis may be easily generalized to more general compact groups (at least to 
SU(n)). 

First of all, we will show that M^ 1 '^ satisfies the commutation relation, 

= if a p 7 M^ r \ (A15) 
Consider the integral 

Let us make a change of integration variable g — > gh. Since the measure is right invariant, 
we have 

j dg (D*(gh)Y D%{gh)Dl{gh) = J dg (D^g))* D%)D%(g). (A17) 
For an infinitesimal h, it leads to 

10 



(Tl,r) 



a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


b 


b 


b 


b 


b 





q 





p+q-l 



p+q-1 



FIG. 1: The Young tableaux method of forming the adjoint (octet) representation in the direct 
product 1Z®1Z. After making p + q — 1 singlet combinations, there are two ways of forming an 
adjoint representation. If q = 0, the way shown on the right is impossible. 



where the independence of the CG coefficients has been used. From this, the commutation 
relation (jA15|) directly follows. 

Next, we show that such a matrix that satisfies Eq. (jA15|) is a linear combination of 
and = d^T^T^ . (Do not confuse with the representation matrix D^j(g).) If the 
Dynkin index of 7Z is (0, q) or (p, 0), they are not independent: is proportional to T^. 
The proof goes as follows. It is easy to show that and D 1 ^ satisfies Eq. (jA15|) . The point 
is that the multiplicity m of 8 in TZ cg> 1Z is at most 2 as we shortly show, so that and 
span the complete set. 

Let us consider the CG decomposition of 1Z ® 1Z using Young tableaux (Littlewood's 
method). (See, for example, Ref. jS^|.) Suppose the Dynkin index of the representation 1Z 
is (p, q), and, without loss of generality, p > q. There are (p + 2q) + (q + 2p) = 3(p + q) boxes 
in the product 7Z<^)71. As is shown in Fig. Q there are at most two ways to form an adjoint 
representation. This explicitly shows that the multiplicity m of 8 in 1Z <S> TZ is at most 2. 

Similarly, we can show that for SU(n) there are at most n — 1 ways of forming an adjoint 
representation from the direct product 1Z®1Z. This number is just the rank of the group, 
rank(G), and it is also the number of invariant tensors. The matrix M^ 1 '^ may be written 
clS 8b linear combination of rank(G) quantities, 



rp(JZ,s) _ rpTZrp'R. rpl 



i7? 



rank(G)) 



(A19) 
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where g a a 1 a 2 ---a s is a real symmetric invariant tensor of SU(n), and therefore T a n ' s ^ is Her- 
mitian. For SU(3), there are two symmtric invariant tensors, 5 a p and <i Q( g 7 . Independence 
of may be examined by defining a rank(G) x rank(G) matrix C st , 

C st id n = T^ s) T^\ (A20) 

Note that the right hand side commutes with for all (3, so that it is proportional to idn 
by Schur's lemma. In general, m = rank(C) < rank(G). We arrange that the first m Ta' s ^s 
are independent. 

We may now write the matrix Ma as a linear combination of independent s, 

m 

MS l ^ = Yv n TS t ^. (A21) 

8=1 

Note that V rs is a regular m x m matrix, but it is not orthogonal. 

The normalization of Ma is determined by the orthogonality condition of the CG 
coefficients, 

{{M^) ] M^) ^ = 5 l3 5 rs . (A22) 
By substituting (|A21jl into the above expression, we get 

{VCV*) n = 6 rs , (A23) 

thus, 

We can write down (|A16J) in terms of T^ ,r \ 

r=l 
s,t=l 

By substituting (|A24|) into this, we finally get the expression, 

/ dg mg))* D%)Dl{g) = ± (C~X (T^) kl (T^) jV (A26) 

For SU(n), £> (8 %) should be replaced by D Ad (g). 
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3. Formulae for the matrix elements diagonal in representation 



Let us now calculate the matrix elements of various operators, which are diagonal in 
representation, by using the formulae derived in the previous subsection. 

For SU(3), there are two symmetric invariant tensors 5 a p and <i Q/ g 7 , we have 

T^) = T?, = d^lfl? = D* (A27) 

The matrix C rs may be written as 

l C 2 C 3 \ 

C=\ 3 , (A28) 

\C 3 D 2 J 

where C 2 and C3 are quadratic and cubic Casimirs respectively, 

C2 = \ (p 2 + q 2 +pq + 3(p + q)) , (A29) 
C 3 = ^(P " l) ( 2 P 2 + 2 1 2 + 5 ^ + 9(p + g + 1)) , (A30) 
for representation 1Z = (p, q), while D 2 = D^D^ may be written as 

D 2 = (~C 2 + ±) C 2 . (A31) 
By looking at the determinant of C, 

det C = -^pq(p + 2)(q + 2)(p 2 + q 2 + 2pg + 4(p + g) + 3), (A32) 
we see that and Z)^ are not independent for p = or q = 0. Actually, 

^ = (A33) 
In this case, the integral (jA26j) gets simpler, 

/ do (omr D { %)Di (g ) = j^^ CTm^v ( a34 > 

When det C^O, and are independent and the matrix C has the inverse, 

_! 1 / C 2 (3 + 4C 2 ) -12C 3 . 

C = Cf(3 + 4C 2 ) - 12C| I _ 12C3 12C2 
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a. Hamiltonian operators 



Here we explicitly use Yr = 1 and present the results for the various flavor SU(3) oper- 
ators appeared in the Hamiltonian. They are expressed in terms of flavor and D™, i.e., 
in the Gell-Mann-Okuboj^l form. 

Let us first show explicitly the simplest case. 



^F 2 S 2 



D 



(8) 
a/3 



*l Sl ) = dhnTZP*(S 2 )P(S 1 ) (| dg (Dl §2 (g)y D^D^g) 



Since 



nS 2 )P( Sl ) £ (c-i) st (T^) F2Fi (T^) 

s,t=l 



l J s ) Sl s 2 ~ 2 5lS2 ' 

WW, = x ( J2 " 5^ - 1) 



(A36) 



(A37) 
(A38) 



we have 



(8) 



(Ds) 



F 2 F! 



C91 + Con 1 ( J 2 — —C\ 



S §1 s 2 , (A39) 



3 4 

where J 2 stands for the spin SU (2) quadratic Casimir. We dropped all the TZ dependence 
on the right hand side for notational simplicity. 

Other matrix elements may be calculated in a similar way. The results are summarized 
as follows, 



+ W F2Fl (c^oP + c^o\ 2) ) ( A4 °) 



S1S2 



14 



where 



O^ 1 

^ X 

o{ 2 



o^ 1 

y 



(2 

^ z 

^ Z 
W 

(2 



2 

2 J ' 



3J' 



Oi-jlJ 



3 1 „ 2 3 2 
8 + 2 C 2 + 3 Cs " 2 J ' 



16 4 2 

6 
3 



3J 4 
1 



Co-- 



C 2 J 2 -(J 2 ) 2 



Co 



1 

3 



-C^J 2 



2\2 



(A41) 

(A42) 

(A43) 
(A44) 

(A45) 

(A46) 

(A47) 
(A48) 



32 8 3 v ' 3 4 6 
The usefulness of these formulae rests on the fact that they are easily calculated for arbitrary 
representation 1Z. 



b. Decay operators 

The matrix elements of the Go operator may be calculated by using Eq. f)A36|) . 

s,t=l 

= p*(s 2 )p(5 1 ){(t q ) F2Fi (cn 1 (7%s 2 + cr, 1 (a)^) 

+ {D a ) F2Fi [C^ (T t ) 5l4 + C^ 1 (Di) §i§2 ) }. (A49) 

Since Tj is the usual spin SU(2) generator, the matrix elements between the states with 
different spins vanish and only the Di terms contribute. When the both states have the 
same spin, from the transformation property, Di is proportional to Tj. By using 




(A50) 



we have for the same spin states (with Yr = 1), 

W §1 s a =A(T i ) §lSa , (A51) 
15 



^F 2 S 2 



D 



(8) 



with 



(A52) 



The matrix elements, diagonal in representation, of the decay operators which are linear 
in F a may be calculated in a similar way. For the G 2 operator, we have 

2 



^F 2 S 2 



D {8) F 



= ns,)P(S l ) £ (C-)„ (7f)„v(*f%* (-7f ^ 



:,t=l 



(A53) 



By using Eqs. (|A37|) and (|A38J) . this can be rewrite as 

*l Sl ) = -^P*(S 2 )P( Sl ){(T a ) F2Fi (ctf + Ctf (j 2 ~\c 2 -\ 



^F 2 S 2 



n( 8 ) p. 



Similarly, we can calculate the matrix elements for the G\ operator, 



(A54) 



K 

F 2 S 2 



(A55) 



s,t=l 



Note that 



di^TpT^ = A —TiT s , 



\C 2 + ~) T, - -L (AT 8 + A^) 



(A56) 
(A57) 



we have 



-P*(5 2 )P(5!)| (T a ) F2Fi (W(A - Tj) + C^ 1 (Qc 2 - ij 2 + |) T 4 - ^ ^ s 

+ (a*) F2Fi (c^a - Tj) + c^ 1 (Qc 2 - i J 2 + ~) ^ - 1a) ) . , }■ (A58) 



When det C = 0, the corresponding formulae become much simpler. 
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APPENDIX B: TABLES OF VARIOUS MATRIX ELEMENTS 



Those matrix elements of the symmetry breaking operators which are needed to calculated 
the baryon masses are given in Sec. 1111 Bl In this Appendix, we summarize other matrix 
elements which are necessary to calculate the mixings to second order. 

These matrix elements are calculated with the method explained in the previous section. 

To second order, the (mainly) octet states can mix with 10, 27 d, 35^, and Q4 d . 
The (mainly) anti-decuplet states can mix with 81^ in addition to them. We there- 
fore need (27 d \Oi\ 27 d ), <35 d | C»< | 35 d > , (35,* |£>i| 27 d ), (64 d \Oi\ 27 d ), (64 d \Oi\ 35 d ), and 
(8T d \Oi\35a). 



TABLE XVIII: (27 d \Oi\ 27 d ) 





7 


X 


y 


z 


w 


(i+i) 


137 


3 137 




137 


2803 


560 


4 560 


1120 


112 


1120^ 


(1,0) 


13 
280 


3 13 

4 280 


29^ 
560 


13 
56 


407 
560^ 


(i-i) 


2 

35 


3 2 

4 35 


13^ 

35 


2 
7 


311 
140 


(0,0) 


13 
70 


3 13 

4 70 


29y/~3 
140 


13 
14 


407 
140 


(1,-1) 


17 
112 


3 17 

4 112 


31^ 
224 


85 
112 


235 
224^ 


TABLE XIX: (35 d \0-\ 35 d ) 


V,Y) 


7 


X 


y 


Z 


W 


(0,+2) 


i 

4 


3 1 
1 i 


8 


9 
1 


8 


(i+i) 


3 

16 


3 3 

4 16 


llv^ 
32 


27 
16 


43^ 
32 


(1,0) 


1 

8 


3 1 
1 S 


16 


9 
8 


17v/3 
16 


(1,-1) 


1 

16 


3 1 

4 16 


7^ 
32 


9 

16 


25^ 
32 


The 


(mainly) decuplet states can mix 


with 27g, 35, 35 f/ , 


64 g , and 81. 


Thus we need 


<27 g |0 4 |27 9 ), (3510*135), 


<35|0 4 |27 g ), 


(35,10*127,), (64 t 


,\Oi\27 q ), (64 


q \Oi\ 35), and 



(8110,1 35). 

Because these matrix elements listed above only contribute to the second order calcula- 



tions, those of Vl and V2 are not necessary, because these operators themselves are of sec- 
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TABLE XX: (35, \0 { \ 27,) 





7 


X 


y 


z 


w 




16\/7 


3 VIE 

4 16-v/7 


32^7 


VToE 

16 


11VE 
32V7 


(1,0) 


8v/7 


3 VE 

4 8^ 


nVE 

16V21 


8 


11VE 
16V2I 


(1,-1) 


16-/7 


3 VIE 

4 16-/7 


nVE 

32V7 


VT05 
16 


11VE 
32V7 


TABLE XXI: (64, \Oi\ 27 ,) 


(I Y) 


7 




y 






v 2 ' 1 / 


5V3 
28 


3 5V3 

4 28 


_5_ 

14 


85^ 
56 


5 

56 


(1,0) 


5^ 
28^ 


3 5^ 

4 28^ 


5VE 

14^6 


85^ 
56^ 


5VE 

56^ 


(1,-1) 


5^ 
28^ 


3 5^ 

4 28\/2 


5 

14^ 


85^ 
56^ 


5 

56^ 


(i-i) 


5\/3 
28 


3 5^ 

4 28 


5 

14 


85V 7 ! 
56 


5 

56 


(0,0) 


3^5 
Uy/2 


3 3^ 

4 14%/2 


VIE 


51^ 
28^ 


VIE 
28^ 






TABLE XXII: ( 


64,10,135,) 








7 


X 


V 

•J 


Z 


W 


(|,+1) 


1 

4^ 


3 1 

4 4^ 


\[i 


3V1 

sVE 


13^ 

sVsE 


(1,0) 


1 

4 V / Ti 


3 1 

4 4VT4 


IT 
V 14 


3^/7 
873 


13-y/Jj 


(1,-1) 


3 


3 3 

4 4^ 


3^ 




39\/3 






TABLE XXIII: 


(81 d \Oi\M d ) 






fry) 


7 


X 


y 


z 


w 


(0,+2) 


2^ 


3 y/3 

4 2^ 


1 


nVs 


1 


(i+i) 


2 

v 7 ^ 


3 2 
4^ 


1 

VwE 


26 
v 7 ^ 


1 

VwE 


(1,0) 


2V7 


3 v 7 ! 

4 2^ 


i_ 

4V7 


13^ 
2^ 


1 

4^ 


(1,-1) 




3 AT 

4 V 35 


1 

2^ 


13^ 
V3~E 


1 

2^ 


ond order. There are, however, extra matrix elements that we 


need to calculate; (35, 


lais), 


(64, |0 


v \8), (64 d |a|10>, 


(81,1a! 10), <35 g 


|a|10), (64 


q \O v \ 10), and (81 \O v 


10). 
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TABLE XXIV: (27 , \Oi\ 27 q ) 





7 


X 


y 


Z 


It' 


(i.+i) 


13 
112 


15 13 
4 112 


224 


13 
56 


11 

32^ 


(1,0) 


1 

56 


15 1 
4 56 


25V 7 ! 
112 


1 

28 


1 


(i-l) 


17 
112 


15 17 
4 112 


145^ 
224 


17 
56 


7 

32^ 


TABLE XXV: (35 1 0*1 35} 


fry) 


7 


;r 


y 


Z 


W 


(i.+i) 


3 

16 


15 3 
4 16 


32 


9 

8 


nV3 

32 


(1,0) 


1 

8 


15 1 

4 8 


3^ 
16 


3 
4 


ll-v/3 
16 


(i-i) 


1 

16 


15 1 
4 16 


23^ 
32 


3 
8 


27^ 
32 


(0,-2) 








5\/3 
4 





-Vs 


TABLE XXVI: (35 |0<| 27 g ) 




7 


;r 


y 


z 


w 


(f.+l) 


VIE 

16%/7 


15 %/l5 
4 16\/7 


32^ 


Vl5 


nVE 
32V7 


(1,0) 


8V7 


15 

4 8^ 




V5 
2V7 


nVE 
16V21 


(1,-1) 


\/l5 
16-/7 


15 v 7 ^ 
4 16^7 


V5 
32^7 


VIE 


17V5 
32^ 






TABLE XXVII: 


(35, | Oil 27, > 






V,Y) 


7 


X 


y 


z 


W 


(1,+1) 


\/5 


15 FT 
4 V 35 


2V7 


V35 


V5 
2V7 


(1,0) 


i 

2V7 


15 1 
4 2V7 


5 

4 V / 2T 


2 

V7 


5 

4v / 2T 



APPENDIX C: "TRADITIONAL" APPROACH TO THE SKYRME MODEL 



In this section, we consider the "traditional" approach to the action ()2.9|) . Namely, we 
calculate the profile function F(r) of the soliton, then all of the Skyrme model parameters 
are determined by the parameters and the integrations involving the profile function. 

It has been known that, in the conventional Skyrme model where 1 is the only symmetry 
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TABLE XXVIII: (64 q \Oi\ 27 q ) 





7 




:r 


y 


z 


W 


(i.+i) 


5-^3 
56 




15 5^ 
4 56 


5 
7 


55^ 
112 


5 

16 


(1,0) 


56 




15 5VE 
4 56 


5^ 
7V3 


55^ 
112 


Sv 7 ! 

16-^/3 


(i-1) 


5V 7 ! 
28V2 




15 5V3 
4 28^ 


5V2 
V7 


55V 7 ! 
56^ 


5 

8^ 


(0,0) 


3VE 
28 




15 3^ 
4 28 


2y/l5 
7 


33^ 
56 


8 


TABLE XXIX: (64 q \Oi\ 35) 




7 




X 


y 




it; 


(1,+1) 


9 

8^ 




15 9 
4 8^35 


3VTE 

4V7 


27^ 
16-v/7 


5l-v/3 
16-v/35 


(1,0) 


3 

8^7 




15 3 
4 8^ 


5V3 
4V7 


45 
16^7 


17V3 
16%/7 


(i-i) 


3 

4^70 




15 3 
4 4^ 


2 V / 14 


9^ 
8^ 


17%/3 
8^ 


TABLE XXX: (81 |£>i| 35) 




7 




X 


y 


Z 


»' 


(f.+l) 


3 

2V35 




15 3 
4 2^ 


4^7 


3^ 
V7 


3^ 
4^ 


(1,0) 


3 




15 3 
4 4^ 


15V3 
8^ 


15 
2^ 


3^ 
8^ 


(1,-1) 


yi 




15 AT 

4 V 35 


3^ 
2s/7 


2^ 
V7 


3 

2^ 


(0,-2) 


3 




15 3 
4 4^ 


15V3 
8^ 


15 
2^ 


3^ 






TABLE XXXI: ( 


35 d |a|8> & (64 d |a 


8) 










(35 d |a|8> 




(64 d |a|8) 




(I,Y) 








^2 


VI 


U 2 


&1) 


3 




1 


3 

28^ 


1 

28^ 


(1,0) 




l 

2V35 




l 

6^ 


V5 
28V2 


75 
84^ 


(|,-1) 












3 

28^ 


1 

28^ 


(0,0) 












9 

28^ 


3 

28^ 



breaking interaction, the physical values of the parameters do not reproduce the baryon 
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TABLE XXXII: <64 d \O v \ 10> & (81 d \O v \ 10> 







/r4j \ci 1 1 n\ 




/si ,,\o I i n\ 

\ OJ -a | *-sv \ - Lu / 


(I,Y) 






V\ 




(0,2) 








56 


56 


(ii) 


3 


1 


9 


3 


28vT0 


28^ 


28^ 


28^ 


(1,0) 


1 


1 

42^ 


3^ 
56 


V3 
56 


(1,-1) 


3 

56 


1 

56 


3^ 
56^ 


V3 
56^ 



TABLE XXXIII: (35, 10), (64 g \O v \ 10) & (81 |O w | 10) 







(35ja| 10) 


(64 g | 


O v \ 10) 


(si i a 


10) 


(I,Y) 




V2 




v 2 


Vl 




(1,1) 




V5 


3^ 




3V5 


Vf> 


6%/2 


18V2 


56^ 


56V2 


112-v/2 


112V2 


(1,0) 




V5 
18-y/2 


28 


VE 

84 


3^ 
112 


112 


(i-i) 








3 

56 


1 

56 




V3 


(0,-2) 














3V5 
112 


75 
112 



mass spectrum. In particular, the best fit value of the pion decay constant F n is typically 
less than one third of the physical value (e.g., 46 MeV[30] for the physical value 184.8 MeV), 
while the Kaon mass (if it is treated as a parameter) becomes quite large (around 800 MeV). 

From the effective theory point of view, this discrepancy can be understood easily. There 
are an infinite number of operators in the action and they contribute to the Skyrme 

model parameters. The conventional Skyrme model ignores all of such contributions from 
the higher order terms. Furthermore, we do not know those coupling constants at all. In 
this paper, we therefore give up to "calculate" the Skyrme model parameters from the 
parameters, and fit them directly to the experimental values. 

An important question is whether the best fit values of the parameters "improve" 
as we systematically take the higher order contributions into account. In this Appendix, we 
address this question. Starting with the action (|2.9[) . we calculate the profile function and 
the Skyrme model parameters. By fitting them to the experimental values of the baryon 
masses, we obtain the best fit values of the %PT parameters. 
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1. Profile function F(r) 



First of all, it is useful to define the subtracted action S su b[U] as 

S sub [U] = S[U] - J rf 4 xTr (M* + M) - L S B 2 J d"xTi (M ] M ] + MM) , (CI) 

so that S su b[U = 1] = 0. 

The profile function F(r) minimizes M[F] defined by 

S sub [U c ] = —M[F] J dt, (C2) 

subject to the boundary conditions, 

F(0) = vr, F(oo) = 0. (C3) 

The minimum is called M c i, 

M d = mm M\F]. (C4) 

F(r) 

The profile function satisfies the Euler-Lagrange equation, 

d dM dM 0. (C5) 



dr d(dF/dr) dF 
By substituting the solution F(r), we obtain M c i, 
4vrF 



M d = / r z dr 

e Jo 



1 |\ ,, 2 2sin 2 Fl sin 2 F f ... 2 sin 2 F 



m5n , , AmBnLz. { , .. 2 2sin 2 F 

+ (1 - cosF) + — ^cosF |(F') 2 + — ^ 



, 8m 2 F 2 L 8 . 2 

+ F 4 2 Sm F 



(C6) 



where we have introduced f = Foer and F' = dF/dr. Note that there are higher order 

contributions and they affect the behavior of the profile function, and hence the values of 

the other parameters. The existence of the higher order interactions affects the values of 

the parameters in two ways; through the behavior of the profile function and the new terms 

from the higher order interactions. 

For a large value of r, F(r) and its derivative behave like F ~ 0, F 1 ~ 0, and thus 

sin F ~ F, cos F ~ 1, we have 

dM 4ttF fl , 8mB L 5 \ ~ 2 w ™ 

1 r F , (C7) 



<9F' e V 4 F o 
dM 4vrF / mF 16m 2 F 2 L 8 x 

W~—\2F^ + W 2 ' ( 8) 
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By solving the Euler-Lagrange equation, the asymptotic behavior, 



F(f) 



-fir 



~ 5 

r 



(C9) 



is obtained, where /i is given by 



F e \ 
F e 



-i , 32mB L8 
1 i pi — 

1 H — - 

16mB 



F 2 



(L 8 - U 



(CIO) 



Note that \i is close to the pion mass but not exactly the same. 



2. Inertia tensor I a p(A) 



After substituting U = AU C A^ into S su b[U], the inertia tensor can be easily read off as 
the coefficients of the terms quadratic in the "angular velocity" u a , 

J d 3 a^Tr [KU%U C + - {X a , Xp}] 

f d 3 xTr M\ a U c - ^UtdiUMXpUc- XpMdiUc]) 
IZoe J 

2L 5 B 2m + m s J ^ ^ + ^ ^ p) 1^ + _ A/?} ] 



F 2 
zJL 

16 



3 3 
L5B0 2m + m s 
3 



rf 3 x(l-cosF) -Tr 



X s {X a Ut\pU c + X p UtX a U c 



+ U\X a U c X p + U\X p U c X a - 2{A a , Xp}) 



L 5 B m s ^_m j ^ 1 ^ ^t AgA) ( [Aq) j/J j/t Uc] + Uc] [/t [Aqj ^ 



+ [A Q , U\] U c [Xp, U\] + [A* U\] U c [X a , U\] )] 
The A-independent part ]L may be written as Eq. (|2.24|) with 



2a 



Ji = — 5— / r df sin F 



71 



h = _ " / r rfr(l-cosF) 



1 + 4 j(F 
1 + (F') 



,, 2 sin 2 F] 327715, 



cos F 



2 , 2sin 2 F , SmB 

L 5 (1 + cosF) 



Ft 



where 



2m + m. 



(Cll) 

(C12) 
(C13) 

(C14) 



53 



has been introduced. 

The ^4-dependent part I' a p(A) has the complicated structure as Eq. ()2.25j) . The parame- 
ters are now expressed as 



647rL B Sn 



CO 



x = — - — ^ 3 ° 5m / r 2 dr cos F sin 2 F, (C15) 
9e -fo jo 



32ttL 5 B 



y = — -8m / r 2 dr sin 2 F, (C16) 

3V3e*F 3 Jo K J 

on t JD POO 

z = "° 5m / f 2 ^(l-cosF)(2-cosF), (C17) 



9e 3 F 3 Jo 
16kL r B 



W = ^Tvi 5m I r 2 dr{l-cosF){2-cosF). (C18) 
S'v 3c Fn 



3. Potential 

The "angular velocity" independent part (except for M c i) of the Lagrangian is the po- 
tential part Eq. (|2.27j) . —V(A). The parameter 7 in the first-order term ()2.28|) is given by 



with 



x ^B 
70 = 5m ^F n 



7 = 7o + 5i7 + $j7> (C19) 



roc 

° / f 2 df(l -cosF), (C20) 
Jo 



327Tj? L 5 - / 2 2sin 2 F . 

,3 l7 = dm - - — / r dr cos F [(F) H — — , (C21) 

3eF J V H 



527 = 5m ^|o^s r f 2 df(1 _ cos2F) (C22) 



The parameter t> in the second-order term ()2.29|) is given by 

1671-56% 



„ = (£ m ) 2 Q 3 8 / f 2 ^(l-cosF)(l-2cosF). (C23) 



4. Numerical calculations 



Starting with the xPT parameters Fq, Bq, e, L5, L§, and quark masses, m and m s , we 
can first calculate F(r) and then, by using it, the Skyrme model parameters, M c i, Ji, 7 2 , 7, 
x, ?/, 2, w, and v. Once these parameters are determined, the baryon masses can be easily 
calculated in the second-order perturbation theory. In order to best fit the xPT parameters, 
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we need to solve reversely. The procedure is similar to that discussed in Sec. IIVI but a bit 
more complicated. In order to simplify the calculation, we make the following things; 



1. The quark masses are fixed. As a reference, we adopt the following values 



m = 6 MeV, m, = 150 MeV. 



(C24) 



Actually, it only fixes the ratio, because the change of the magnitude can be absorbed 
in Bn. (Note that M. appears only in the combination BnAi.) The ratio is better 
determined experimentally and known to bep, 166] 

2(m s /m d ) 



(m u + m d )/2 l + (m u /m d ) 
which is close to the value 150/6 = 25. 



25.8 



(C25) 



2. The values of L5 and L$ are fixed. When we vary these parameters too, we find that 
numerical calculation becomes very unstable. In reality, all of the formulation in this 
paper assumes these parameters to be small. In searching the "valley" numerically, this 
assumption is often ignored, and we believe that it is the reason of the instability. In- 



stead, we fix these parameters to be the central values determined experimentally 



4l|, 



L 5 = 1.4 x 10" 3 , L s = 0.9 x 10^ 3 . (C26) 

There is another important point. As Yabu and AndoQ] discussed, there is a kind 
of "zero-point energy" contribution universal to all of the calculated baryon masses. This 
contribution may be calculated as the symmetry breaking effects to the fictitious (unphysical) 
singlet baryon mass, 



M„ 



7 



1 ) + V[ 1 



M 8 - Mi 
7 ,v _ 7 2 / 2 
2 2 48 ' 

where we have introduced 

M 1 = M a + - 

o 

We subtract it from all of the calculated masses 



1 2 

h~T 2 



(C27) 



(C28) 
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Note that, in the effective theory approach, this contribution is renormalized in the 
parameter Ma, and, therefore, does not need to be considered separately. It has been 
implicitly taken into account. 

Our numerical results are 

F = 82.7 MeV, B = 2697, e = 4.51, (C29) 
which lead to the following values for the physical parameters, 

F n = 91.4 MeV, M n = 185.3 MeV, M K = 867.2 MeV (C30) 
and the baryon masses for these values are given in Table IXXXIVI These values should be 

TABLE XXXIV: Baryon masses for the best fit values in the "traditional" approach. The results in 
the first low is with the higher order contributions, and in the second low without. The "zero-point 
energy" contribution has been subtracted. 



Baryon 


M N 


M s 


M s 


M A 


M A 


Ms* 


Ms* 


Ma 


M 


M^ 


Mass™ (MeV) 


915 


1287 


1411 


1116 


1185 


1358 


1518 


1666 


1563 


1965 


Mass„, /o (MeV) 


898 


1269 


1405 


1116 


1118 


1321 


1501 


1660 


1639 


1948 



compared with those calculated with L5 = L§ = 0, that is, those without the contributions 
from higher order terms. 

F = F n = 58.2 MeV, B = 7825, e = 4.06, (C31) 

which lead to 

M w = 306.4 MeV, M K = 1104.9 MeV. (C32) 

The baryon masses are also given in Table IXXXIVl 

It is interesting to note that the values of the physical parameters shift in the right 
direction. Even though these values are still far from the experimental values, we think that 
this is an explicit demonstration that our basic strategy is right. 

APPENDIX D: SYMMETRY BREAKING INTERACTIONS IN THE CHIRAL 
QUARK-SOLITON MODEL 

In this Appendix, we perform a similar "best fit" analysis with the symmetry breaking 
terms (jl.4|) which appear in the xQSM. In the derivation of these terms |67(, they are related 
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to the 7r-N a term, soliton moments of inertia, and so on, but we ignore this fact and just 
treat the couplings as free parameters. The reason is that it makes the comparison with our 
approach transparent and reveals how the xQSM predictions depend on the detailed form 
of the parameters. 

1. Best fit to the baryon masses 

To obtain the masses in the second order perturbation theory, we need the matrix elements 
of H® pp . The matrix elements of are given in Sec. IIII Bl and in Appendix El an d those 
for Y are trivial. It is only Oppp = D Si F l of which the matrix elements need to be 
calculated. We first present the matrix elements for the spin J — \ states in Tables IXXXVI 



and lXXXVII 

TABLE XXXV: The matrix elements of Odpp, (R-i \Odpp\ for the spin doublet states, which 
are abbreviated as (1Zi,lZj). 





(8,8) 


(8,10) 


(8,27 d ) 


(10,10) 


(10,27,,) 


(10,35*) 


(0,+2) 











8 







(i+i) 


20 




l 


^3 
16 


7 

16^ 


3^ 


(1,0) 


3^ 
20 


V3 


1 

lOv 7 ^ 





7 

sVIE 


V3 


(i-i) 


5 





1 

lOv^ 






















V3 
16 


7 

16a/2 


VIE 

16\/l4 


(0,0) 


373 
20 





20 












Those for the spin J = § states in Tables IXlOCvTTl and IXlQCVTHl 

Most of them have never been given in the literature. (Incidentally, the matrix elements 
diagonal in representation are given by similar formulae to (|A40J) given in Appendix 1X1 with 

= -J\ (Dl) 
0%„ = - (ic s + \C 2 + ij» + 1) . (D2) 

See Sec. IA 31 for the notation.) 

By using these matrix elements, we can calculate the baryon masses, and by best fitting 
the calculated values to the observed ones, we can determined the parameters a, (3, and 7 
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TABLE XXXVI: The same as in Table IXXXVI 



(I,Y) 


(27 d , 27 d ) 


(27 d ,35 d ) 


(27 d ,64 d ) 


(35 d ,35 d ) 


(35 d ,64 d ) 


(35 d ,81 d ) 


(0,+2) 











8 





3y/3 
ov ' 


(i+i) 


71 V3 
1120 


5^ 


5 

56 


3^ 
32 


3^ 

ov OO 


3 


(1,0) 


19>/3 
560 


5^ 
16^/21 


5\/5 
56^ 


16 


8VT4 


_3a/3 
8%/7 


(i-i) 


ll-y/3 
70 





5 

56 











(l,-i) 


224 


5^ 
32^ 


5 

56^ 


32 


9^ 
8^ 


3^ 
4^ 


(0,0) 


19-^3 
140 





VIE 
28\/2 












TABLE XXXVII: The matrix elements of Odpp, (T^-i \Odpp\ T^j), for the spin quartet states, which 



are abbreviated as (TZi,TZj). 





(10,10) 


(10,27,) 


(10,35) 


(27 ff ,27 4 ) 


(27„35) 


(i+i) 


5^ 
16 


5^ 
16-^/2 


5^ 
16^14 


65^ 
224 


32 


(1,0) 





5 


_VTS 
8^ 


5^ 
112 




(i-i) 


5yJ 
16 


5 

16-v/2 




85^ 
224 


32 


(0,-2) 


5%/3 
8 





VIE 

8y/7 









TABLE XXXVIII: The same as in Table IX XX VI II 





(27 9 ,35 9 ) 


(27 5 ,64 9 ) 


(35,35) 


(35,64 g ) 


(35,81) 


(1,+1) 


VE 

2v/7 


5 

112 


32 


3^ 
16^ 


3^ 
4^ 


(1,0) 


5 


5^ 
112-v/3 


3^ 
16 


V21 
16 


3V3 

8V7 


(i-i) 





5 

56^ 


ll\/3 
32 


8-y/IO 


3 

2^ 


(0,-2) 








2 





3V3 
8V7 



as well as M c /, Jj, and J2. The procedure is the same as that employed in Sec IIVI so that we 
do not explain it again. The best fit set of parameters is 

M d = 837 MeV, If 1 = 163 MeV, I^ 1 = 394 MeV, 
a = -554 MeV, (3 = -40.9 MeV, 7 = 42.0 MeV, (D3) 

which leads to the masses given in Table IXXXIX| with x 2 — 6.4 x 10 1 . 
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TABLE XXXIX: Baryon masses calculated by using the best fit set of parameters (|D3|) with the 
xQSM Hamiltonian. 



(MeV) 


N 






A 


A 


S* 




n 


e 


4> 


Mf pp 


942 


1206 


1335 


1116 


1227 


1383 


1531 


1672 


1538 


1868 



Considering the number of parameters is small, the fit is very good. It is also remarkable 
that the parameters have expected magnitudes. The best fit values predict the masses of 
the other members of pentaquarks, 

M w = 1668 MeV, M E , = 1777 MeV. (D4) 

The results are not very different from the ones in ()4.4I) . 

2. Decays 

Next we turn to the calculation of the decay widths. Again the procedure is the same as 
in Sec [V] Actually the necessary matrix elements are the same. The only difference comes 
from the mixings. The mixing coefficients which correspond to our results in Tables IXTllXIIl 



and ETUI are given in Tables KlJ EUl aud lxTTTl 

TABLE XL: Mixing coefficients for the (mainly) octet states with the xQSM symmetry breaking 
terms. They correspond to those in Table IXTl 



Hi 


N 


s 




A 


8 


1 (-0.030) 


1 (-0.027) 


1 (-0.010) 


1 (-0.015) 


To 


0.202 (-0.023) 


0.202 (0.030) 


0(0) 


0(0) 




0.140 (0.003) 


0.114 (0.022) 


0.140 (0.020) 


0.171 (0.008) 




(0.022) 


(0.021) 


0(0) 


0(0) 


64 d 


(0.010) 


(0.008) 


(0.010) 


(0.014) 



One can easily see that the mixings are much smaller for the xQSM breaking terms. 
From the goodness of the fit to the masses and the smallness of the mixings, one may think 
that the perturbative treatment of the symmetry breaking terms is justified for this model. 
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TABLE XLI: Mixing coefficients for the (mainly) decuplet states with the xQSM symmetry 
breaking terms. They correspond to those in Table IXI*T1 



TZi 


A 


s* 




Q 


10 


1 (-0.105) 


1 (-0.060) 


1 (-0.026) 


1 (-0.005) 


27 i 


0.451 (0.005) 


0.329 (0.007) 


0.202 (0.006) 


0(0) 


35 


0.081 (0.022) 


0.103 (0.023) 


0.109 (0.019) 


0.103 (0.013) 


35, 


(0.058) 


(0.027) 


0(0) 


0(0) 


64, 


(0.026) 


(0.025) 


(0.017) 


0(0) 


81 


(0.004) 


(0.006) 


(0.007) 


(0.006) 


TABLE XLII: Mixing coefficients 


for the (mainly) anti-decuplet states with the xQSM symmetry 


breaking 


r terms. Thev correspond to those in Table IXIIII 






TU 


e 


N' 


£' 


4> 


8 


0(0) 


-0.202 (-0.009) 


-0.202 (-0.049) 


0(0) 


TO 


1 (-0.009) 


1 (-0.036) 


1 (-0.044) 


1 (-0.033) 


27 d 


0(0) 


0.105 (-0.037) 


0.171 (-0.029) 


0.234 (0.005) 


35 d 


0.131 (0) 


0.139 (0.008) 


0.131 (0.016) 


0.104 (0.019) 


64 d 


0(0) 


(0.012) 


(0.018) 


(0.019) 


8T d 


(0.009) 


(0.010) 


(0.009) 


(0.006) 



With these mixing coefficients, the decay widths are readily calculated, by using our width 
formula ()5.16|) . First we present the best fitted values of Go, G\, and G2 in Table IXLIIII 

These values lead to the decay widths for the (mainly) decuplet baryons given in Ta- 
ble KEEI 

Finally, we obtain the various decay widths for the (mainly) anti-decuplet baryons, which 
are summarized in Table IXLVl 

Despite the good perturbative behavior in the baryon mass fitting, the decay widths vary 
considerably. In most cases, the second order results are very different from the first order 
ones. Even though perturbation theory seems to work good, some negative values appear 
when we expand the amplitudes. In particular, the width of 6 + is predicted much larger 
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TABLE XLIII: Coefficients of the decay operators for the xQSM symmetry breaking terms. The 
notation is the same as in Table IXVl 





(a) 


(a') 


(b) 


(6') 


Go 


5.33 


4.11 


3.38 


3.95 




9.98 


13.7 


21.4 


17.5 


x 2 


43.1 


14.6 


23.9 


36.0 


G 2 


0.08 


0.06 


0.04 


0.05 



TABLE XLIV: Decay widths for the (mainly) decuplet baryons with the coupling constants in 



Table IXLllll 



(MeV) 


K 


r (a ) 


P(a') 


T(6) 


1>) 


A -> Nvr 


1.47 


90.2 


106 


105 


95.8 


£* Avr 


1.18 


34.1 


33.1 


32.9 


33.5 


£* -> Svr 


0.26 


4.69 


5.85 


6.63 


6.19 




0.49 


12.9 


12.7 


13.8 


13.7 



than the experimental values. 
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TABLE XLV: Predictions for the decay widths for the (mainly) anti-decuplet baryons with the 
xQSM symmetry breaking terms. 



(MeV) 


i\. 


1(a) 


F, , 

l(a') 


F, . 
1(6) 


F, , 
1(6') 




-»■ NK 


1 01 


Oo.U 


ZoU 


797 
/ Z / 


/1 90 

4zy 


N' - 


-> N7T 


1 s 9 

lo.Z 


91 n 

Z1U 


ooy 


1 70"? 

1 ( yo 


1 n^n 
lUoU 


TV Tf 

N' - 


-> Nt? 


4.0Z 


-U.oU 


9fi 9 
ZD.Z 


1 n/i 

1U4 


Q 

oo.y 


N' - 


■4 Air 


^ fid 

o.oy 


n 
u 
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